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various  sound  sources  to  manage  noisy  operations  on  a 
daily  basis.  This  study  developed  the  Lookup  computer 
program  and  data  base  to  provide  rapid  estimates  of  out¬ 
door  noise  levels  from  a  variety  of  sound  sources. 

Lookup  accesses  a  data  base  of  archived  results 
(requiring  about  5  MB  disk  space)  from  typical  situations 
rather  than  performing  fresh  calculations  for  each  consul¬ 
tation.  Initial  timing  tests  show  that  Lookup  can  predict 
the  sound  levels  from  a  noise  source  at  distances  up  to 
20  km  in  1  second  on  a  DOS-compatible  personal 
computer  (PC).  This  report  includes  the  Lookup  program 
source  code,  and  describes  the  required  input  for  the 
program,  the  contents  of  the  archival  data  base,  and  the 
program  output.  Lookup  was  written  to  compile  with 
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ESTIMATING  ATTENUATION  AND  PROPAGATION  OF  NOISE  BANDS  FROM  A  DISTANT 
SOURCE  USING  THE  LOOKUP  PROGRAM  AND  DATA  BASE 


I  INTRODUCTION 


Background 

Unavoidable  noise  generated  by  military  training  activities  can  disturb  the  surrounding  community 
and  become  a  source  of  complaint.  Military  planners  must  be  able  to  quickly  and  accurately  predict  noise 
levels  at  distant  points  from  various  sound  sources  to  manage  noisy  operations  on  a  day-to-day  basis.  Predict¬ 
ing  sound  propagation  outdoors  for  broadband  noise  sources  or  for  long  distances  can  require  considerable 
computational  effort.  A  number  of  available  research-grade  computer  codes  can  predict  sound  propagation 
from  outdoor  sources,  all  of  which  involve  some  kind  of  limitation  (in  accuracy  or  model  assumptions)  or 
expense  that  make  them  impractical  for  routine  use.  For  example,  the  Fast  Field  Program  and  the  Parabolic 
Equation  methods  give  very  precise  solutions  to  the  wave  equation,  but  with  a  tremendous  expense  in  central 
processor  unit  (CPU)  time. 

In  fact,  such  expensive  or  inefficient  programs  may  not  be  necessary  for  routine  sound  propagation 
prediction.  Even  with  these  full-wave  acoustical  models,  there  is  still  a  need  to  ensure  the  accuracy  of 
calculated  results  by  comparing  them  with  large  amounts  of  measured  propagation  data.  Moreover,  many 
situations  in  outdoor  propagation  studies  recur.  For  example,  the  micro-climate  near  the  ground  can  show 
diurnal  cycles  whose  sound  speed  profiles  assume  “typical”  shapes  and  thus  produce  similar  propagation;  the 
listeners  are  near  the  ground;  the  sources  are  stationary  or  slow-moving  when  compared  to  the  speed  of 
sound,  and  so  on.  It  may  be  more  practical  to  archive  results  from  typical  situations  and  access  them  through 
computational  lookup  tables,  rather  than  to  specify  new  parameters  and  do  fresh  calculations  for  each 
consultation. 


Objective 

The  objective  of  this  work  was  to  devise  a  method  to  predict  sound  propagation  of  noise  levels 
resulting  from  military  training  activities,  with  the  accuracy  of  full-wave  acoustical  models,  in  near-real  time. 


Approach 

Many  thousands  of  propagation  cases  were  calculated,  and  the  results  were  organized  and  tabulated 
for  easy  access.  The  many  parameters  that  influence  propagation  were  systematically  varied  to  approximately 
cover  most  situations;  i.e.,  the  source  height,  receiver  height,  receiver  range,  frequency,  wind  speed  and 
direction,  temperature  gradient,  and  ground  impedance  were  varied  in  successive  computations  to  build  a  data 
base  of  typical  results. 

Two  computer  programs  were  written  to  access  the  lookup  table  and  produce  estimates  of  the  received 
sound  level.  The  Lookup  program  accepts  an  input  file  that  describes  the  acoustical  source  and  environment 
of  interest,  and  predicts  the  sound  spectrum  and  frequency-weighted  sound  levels  at  selected  locations.  This 
prediction  is  formed  from  the  best-match  case  in  the  lookup  table,  which  is  adjusted  to  account  for  differences 
between  the  stored  and  desired  cases.  For  example,  the  correction  for  attenuation  due  to  molecular  absorption 
is  such  an  adjustment.  Sound  Propagation  (SP)  is  the  second  program,  which  assembles  information  about 
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the  propagation  problem  at  hand,  displays  and  manipulates  it  graphically,  prepares  input  for  Lookup,  and 
retrieves  and  displays  the  predictions. 

Note  that  Lookup  is  a  user-unfriendly  program  that  deals  only  with  the  mechanics  of  the  calculations 
and  the  efficient  access  to  the  tables.  (It  runs  a  typical  case  in  under  1  second  on  the  PC.)  SP  is  the  user- 
friendly  graphical  user  interface  that  makes  it  easy  to  enter,  display,  and  edit  cases  of  interest,  and  quickly 
view  the  results. 


Technical  Requirements 

Lookup  was  written  to  compile  with  MS-Fortran,*  and  will  run  under  DOS  on  any  IBM  compatible 
with  640  kB  random  access  memory  (RAM).  Lookup  also  conforms  to  ANSI  1978  standard  Fortran  and  will 
run  under  the  Unix  operating  system.  The  Lookup  data  base  requires  about  5  MB  of  free  disk  space. 


Mode  of  Technology  Transfer 

The  Lookup  program  and  data  base  will  be  forwarded  to  the  U.S.  Army  Environmental  Hygiene 
Agency  (AEHA)  for  fielding  in  the  “Weather-Based  Noise  Now-Casting  System.” 


*MS-Fortran  is  a  product  of  Microsoft  Corp.,  Redmond,  WA. 
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2  The  Lookup  Program 


Purpose 

The  purpose  of  the  Lookup  program  is  to  provide  a  rapid  estimate  of  sound  levels  at  a  distance 
from  the  sound  source  using  precalculated  tables,  without  incurring  fte  expense  of  a  lengthy  specific 
calculation  for  the  problem  at  hand.  Source  code  for  Lookup  is  included  in  Appendix  A. 


Lookup  Input 

The  Lookup  input  file  contains  all  of  the  physical  parameters  needed  to  specify  an  acoustic 
propagation  problem,  such  as  the  source  and  receiver  heights,  wind  speed  and  direction,  type  of  ground 
surface,  etc.  Also  included  in  the  file  is  the  name  of  the  Lookup  table  data  base  directory  and  the  name 
of  the  index  file  for  the  table.  Appendix  B  shows  a  sample  Lookup  input  file,  which  is  composed  of  a 
strict  format,  for  ease  in  checking  input  errors.  This  input  format  is  fairly  easy  for  the  user  or  the 
computer  to  read,  but  more  difficult  for  the  user  to  write.  Each  major  entry  or  group  of  entnes  is 
preceded  by  a  header  line  describing  the  entry.  These  headers  are  required  by  the  Lookup  program,  and 
most  of  the  entries  are  required,  although  some  are  of  variable  length. 

Version  Number 

The  first  header  and  entry  specify  the  input  file  version  ID  number.  This  number  ensures  that  any 
future  changes  in  the  input  file  foimat  can  be  recognized  by  corresponding  upgrades  of  the  Lookup 
program. 

Directory  Name 

The  next  header  and  entry  pair  contain  the  directory  name  for  the  excess  levels  database  and  index 
file.  The  directory  name  should  be  chosen  to  be  the  greatest  common  path  name  for  the  files  in  the  data 
base.  Note  that  the  data  base  may  extend  over  several  directories,  but  the  directory  name  should  be 
common  to  all  of  the  data  files.  For  example,  as  initially  distributed,  the  Lookup  data  base  has  four 
directories: 

C:\LOOKUPvSRCA 

C;\LOOKUFNSRCB 

C:\LOOKUPvSRCC 

C:\LOOKUIASRCD 

the  Lookup  index  file  name  is: 

C:\LOOKUFMNDEX10.nL 

and  typical  entries  in  the  index  file  are: 

SRCA\ASGP10.AVE  ... 

SRCCXCHGPIO.AVE  ... 

SRCDXDSGPIO.AVE  ... 
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The  directory  name  sent  to  Lookup  should  read: 

C:\LOOKUP\ 

The  directory  name  should  start  on  the  first  character  position  on  the  line  following  the  header. 

Index  File  Name 

The  next  header  and  entry  include  the  Lookup  index  file  name  for  the  excess  levels  table. 
Appendix  C  shows  an  excerpt  from  the  index  file.  The  index  file  contains  a  line-by-line  listing  of  the  file 
names  of  all  of  the  files  in  the  data  base,  paired  with  the  propagation  parameters  for  each  type  of  case 
examined.  The  parameters  include  the  source  and  receiver  heights,  type  of  ground  surface,  etc.  The  name 
of  the  index  file  and  the  file  names  listed  in  the  index  file  should  be  specified  relative  to  the  common 
directory  name,  as  above.  The  index  file  name  should  be  printed  starting  with  the  first  position  on  the 
line  following  the  header.  In  the  example  from  Appendix  A,  the  index  file  name  is  INDEX  lO.FIL  . 

Source  Height 

The  source  height  is  the  height  of  the  source  above  the  ground  surface  (in  meters).  If  the  ground 
is  rugged  or  sloping,  the  source  height  used  is  the  distance  along  a  perpendicular  from  a  local  average 
surface  plane.  The  initial  release  of  Lookup  cannot  form  predictions  for  buried  sources,  and  thus  will  not 
accept  negative  source  heights.  If  the  source  falls  below  the  average  local  surface  height,  but  is  still 
“above”  ground,  the  source  height  should  be  set  equal  to  zero.  Negative  source  heights  will  cause  Lookup 
to  halt. 

Receiver  Height 

All  conditions  that  apply  to  the  source  height,  apply  similarly  to  the  receiver  height. 

Receiver  Range 

The  horizontal  distance  from  source  to  receiver  is  the  “receiver  range”  and  must  be  entered  in 
meters.  If  zero  distance  is  given  as  input,  a  calculation  will  be  made  for  each  range  in  the  data  base,  for 
a  total  of  44  distances.  At  each  distance,  only  the  overall  levels  will  be  reported,  and  no  spectral 
information  will  be  provided.  If  a  positive  distance  is  specified  instead,  a  different  set  of  information  will 
be  given  on  output.  In  this  case,  only  a  calculation  for  that  range  will  be  performed.  Also,  an  estimate 
of  the  receiver  1/3-octave  spectrum  will  be  reported  in  addition  to  the  overall  levels  at  that  distance.  A 
negative  receiver  range  will  cause  Lookup  to  halt.  For  the  format  of  the  output,  see  the  section,  “Lookup 
Output.”  Also  see  “Method  of  Calculation”  for  information  on  the  behavior  of  Lookup  with  regard  to 
receiver  range. 

Receiver  Azimuth 

The  header  and  entry  for  receiver  azimuth  gives  the  direction  of  the  receiver  from  the  sound  source. 
The  units  of  azimuth  are  degrees,  and  are  to  be  measured  from  zero  degrees  north.  Units  of  azimuth 
should  increase  with  clockwise  rotation,  looking  down  on  the  site  from  above.  The  receiver  azimuth  is 
only  used  in  computing  the  wind  component  in  the  direction  of  propagation,  but  must  be  included  in  the 
input  file,  even  if  the  wind  profile  is  unknown  or  unspecified. 
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Ground  Classification 

The  next  entry  is  the  porosity  classification  of  the  surface.  In  the  Lookup  data  base  the  ground  has 
been  classified  on  a  continuous  scale  from  hard  (class  0)  to  porous  (class  1)  as 
ground  between  the  source  and  receiver.  Note  that  not  all  types  of  ground  surfaces  should  be  classified 
fn  this  way  (e  g  ,  snow-covered  ground  should  not),  but  this  scheme  could  be  expanded  for  more  types 
of  surface!  by  Adding  parameters.  Ground  classification  values  outside  of  the  interval  (0,1)  will  cause 

Lookup  to  halt. 

Ground  Roughness  Height 

The  next  entry  is  the  estimate  of  the  roughness  (in  meters)  of  the  ground  surface  (rms).  The  ground 
roughness  is  used  only  in  estimating  the  thickness  of  the  roughness  layer  of  the  wind 
r  rTnimum  roughness  height  for  very  flat  ground  is  about  0.1  m.  Nonpositive  roughness  heights  will 

cause  Lookup  to  halt. 

Temperature  Profile 

The  variation  of  temperature  with  height  is  call  the  temperature  profile,  which 
a  list  under  the  temperature  profile  header.  On  each  line,  the  height  of  the  sensor  should  precede  the 
temperature  measurement.  The  sensor  height  should  be  supplied  in  meters,  and  the  temperature  in  degre 
ceXra^  Speratures  outside  of  the  interval  (-50  T,  50  “C)  will  cause  Lookup  to  halt  A  negahve 
height  and  arbitrary  temperature  should  be  given  to  signal  the  end  of  the  temperature  profile  list.  No 
only  the  first  100  height  and  temperature  pairs  are  used  by  Lookup;  the  rest  are  s^ply  dis  . 
values  for  the  temperature  list  are  substituted  if  there  are  too  few  values  to  construct  a  temperatu 

gradient  over  height. 

Humidity  Profile 

The  next  entry  in  the  input  file  is  the  humidity  profile  as  a  function  of  height.  The  profile  should 
be  supplied  as  in  the  temperature  profile  entry,  with  the  height  first,  and  the  relative  ° 

on  the  same  line.  The  relative  humidity  should  be  entered  in  umts  of  percent.  As  for  the  temperature 
profile,  only  the  first  100  entries  of  the  humidity  profile  will  be  used;  the  rest  are  discarded. 

Values  of  percent  relative  humidity  outside  of  the  interval  (1,100)  will  cause  Lookup  ^ 

molecular  absorption  routine  is  less  accurate  at  very  low  humidities,  and  excluding  the  interva  (  ,  )  e  p 
detect  errors  that  might  occur  when  relative  humidity  has  not  been  properly  converted  to  a  percentage. 

Wind  Profile 

The  wind  profile  is  entered  below  the  humidity  profile  and  follows  a  similar  fonnat.  The  sei^r 
height,  wind  speed,  and  wind  direction  for  each  wind  sensor  are  listed  on  each  line  of  the  profile.  The 
profile  list  should  end  with  an  entry  of  negative  height  and  with  arbitrary  (nonempty)  values  for  wind 
speed  and  direction.  The  wind  speed  and  direction  are  expected  to  be  expressed  in  umts  of  meters/sec 
and  degrees  from  north,  respectively.  The  wind  direction  should  be  entered  in  e  inpu 
meteorological  wind  direction,  i.e.,  the  direction  the  wind  vane  points. 

Source  Level  Spectrum 

The  source  spectrum  should  appear  below  the  meteorological  profiles.  The  source  spectrum  should 
be  entered  as  1/3-octave  band  levels  (dB)  starting  with  band  number  zero  (1  Hz)  and  ending  with  band 
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number  43  (20  kHz).  There  must  appear  values  for  each  band  index.  Unused  bands  should  be  set  equal 
to  -100.0  dB. 

To  obtain  only  the  spectral  decay  from  the  source,  set  all  of  the  bands  equal  to  zero  dB,  the 
reference  distance  (described  below)  to  1  m,  and  specify  the  range  of  interest  (rather  than  letting  range 
equal  zero). 

The  input  spectrum  should  be  flat-weighted,  although  Lookup  will  operate  with  and  perform 
calculations  on  a  frequency-weighted  spectrum.  It  should  be  understood  that  a  frequency-weighted  input 
spectrum  will  not  obtain  proper  values  for  the  A-  and  C-weighted  overall  levels.  The  output  spectra  will 
have  the  same  frequency  weighting  as  that  applied  to  the  input  spectrum. 

The  source  spectmm  should  also  be  taken  from  a  far-field  measurement  and  taken  at  a  distance  at 
which  the  nonlinear  effects  of  propagation  are  weak.  For  example,  measurements  should  be  made  about 
10  wavelengths  from  the  source,  and  at  sound  levels  below  160  dB.  The  effects  of  the  ground  reflection 
are  assumed  to  be  removed  from  the  source  spectrum. 

Reference  Distance 

The  reference  distance  entry  for  the  source  spectmm  follows  the  spectmm  in  the  input  file.  The 
distance  should  correspond  to  the  distance  at  which  the  source  spectmm  was  measured. 

User  Weighting  Spectrum 

The  user  weighting  curve  is  applied  to  the  receiver  spectmm,  and  the  user-weighted  spectrum  is 
summed  by  Lookup  to  give  an  overall  user-weighted  level.  The  user  weighting  spectmm  follows  the  same 
format  as  the  source  spectmm,  with  some  changes.  Only  inactive  (-100.0  dB)  source  bands  can  have 
corresponding  unspecified  bands  in  the  user  curve.  All  other  bands  supplied  from  the  user  curve  are  used. 
There  must  appear  user-curve  values  for  all  1/3-octave  bands  zero  through  43.  The  user  weighting 
spectmm  concludes  the  Lookup  input  file. 


Lookup  Output 

The  output  “file”  is  written  by  Lookup  to  the  standard  output  and  may  be  saved  as  a  file,  viewed 
on  the  display  screen  (as  numbers  and  text)  or  sent  to  another  program  for  graphical  display  or  for  further 
processing. 

The  output  file  will  contain  an  error  report  if  an  error  occurred  during  the  calculation.  If  there  were 
no  errors,  a  report  is  given  that  describes  the  Lookup  parameters  used  for  searching  the  excess  levels  data 
base  and  the  best  match  found  from  the  index.  Appendix  D  shows  the  contents  of  a  sample  output  file. 

The  predictions  for  the  overall  levels  will  be  printed  next.  If  a  particular  range  is  selected  (input 
receiver  range  positive),  the  range,  overall  level,  overall  A-weighted  level,  overall  C-weighted  level,  and 
overall  user-weighted  level  are  printed  on  one  line.  If  all  of  the  ranges  are  selected  (input  receiver  range 
zero),  one  such  prediction  is  made  for  every  range  in  the  excess  levels  data  base  (from  1  m  to  20  km). 
If  a  specific  range  is  selected,  the  predicted  1/3-octave  spectmm  at  the  receiver  will  be  printed. 
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Lookup  Data  Base 

The  excess  levels  data  base  is  composed  of  an  index  file  and  many  files  of  excess  level  spectra  as 
a  function  of  distance. 


Each  excess  spectrum  file  has  a  list  of  range  numbers  and  spectrum  pairs  from  some  starting  range 
(say  1  m  or  10  m)  to  some  final  range  (perhaps  1  km  or  20  km).  The  range  numbers  ^e  defined  as  1^0 
times  the  common  logarithm  of  the  range.  The  range  numbers  advance  by  one  for  each  spectrum  in  the 
list,  such  that  the  corresponding  ranges  advance  through  the  list  by  a  geometnc  progression  in  1/3-octave 

multiples. 


The  spectrum  associated  with  each  range  represents  the  expected  change  in  level  between  the  source 
and  receiver,  due  to  such  effects  as  ground  reflection,  ground  absorption,  refraction,  ^^ff^'^^tion  and  ducted 
propagation  of  sound.  Effects  of  refraction  due  to  changes  in  the  atmosphenc  properties  wi*  heiglrt  are 
LcouLd  for.  The  spectra  do  not  contain  the  effects  of  spherical  spreading  and  absorption,  which  are 

accounted  for  separately  in  the  Lookup  program. 


The  spectra  pertain  to  1/3-octave  bands  of  noise,  rather  than  pure  tones,  although  they  were 
constructed  from  a  pure  tone  propagation  model.  There  is  also  an  implicit  assumption  that  the  source 
spectrum  is  roughly  “flaf  ’  within  any  given  band. 


The  excess  levels  files  were  constructed  from  a  pure-tone,  full-wave,  two-way  wave  equation 
solution,  using  the  Fast  Field  Program  (FFP).  The  FFP  was  set  to  run  a  variety  of  propagahon  problems 
for  all  of  the  1/3-octave  band  center  frequencies.  The  FFP  provides  a  prediction  of  the  excess  level  at  a 
number  of  equally  spaced  distances  between  source  and  receiver. 


The  predictions  from  the  FFP  (or  any  other  wave  model)  usually  vary  in  level  with  distance  due 
to  wave  interference  and  refraction.  (In  other  words,  the  levels  vary  in  response  to  all  of  the  parameters 
we  specify.)  In  some  situations  the  predictions  show  considerable  variation  or  oscillation  with  range, 
circumstances  that  produce  strong  interference  patterns,  it  is  debatable  whether  the  structure  has  any 
practical  significance;  in  a  perfectly  time-invariant  and  precisely-known  environment,  perhaps  it  is 
significant.  In  a  real  environment,  not  one  atmospheric  measurement  can  be  made  with  great  precision, 
nor  will  any  of  the  parameters  remain  invariant  with  time.  Furthermore,  one  seldom  encounters  sources 

that  produce  only  pure  tones. 


As  stated,  the  full  wave  solutions  as  a  function  of  range  do  contain  “structure,”  and  some  of  toe 
rapid  variation  in  range  is  undesirable.  The  more  significant  features  of  toe  predictions  are  those  that 
persist  over  great  distance  or  those  that  do  not  vary  rapidly  with  frequency  or  from  slight  changes  in  toe 

atmosphere. 

Within  any  frequency  band,  toe  band  power  spectrum  is  toe  frequency  integral  of  toe  absolute 
square  of  toe  source  and  transfer  function  product.  If  either  toe  source  spectrum  or  toe  transfer  function 
is  relatively  flat  within  toe  band,  the  integral  can  be  simplified  by  replacing  the  source  spectrum  by  its 
rms  average.  The  integral  can  be  further  simplified  when  toe  oscillating  part  of  toe  transfer  function 
(oscillations  with  respect  to  range)  are  due  to  positional  changes  over  an  interference  pattern.  In  this  case, 
a  similar  oscillating  pattern  can  be  found  by  holding  toe  range  constant  and  sweeping  through  toe 

frequency. 


The  propagation  model  provides  toe  signal  level  at  a  particular  frequency  over  many  ranges,  but 
of  interest  here  is  toe  level  at  a  particular  range,  averaged  over  frequency.  The  average  over  toe  frequency 
interval  should  be  well  approximated  by  toe  average  over  toe  same  number  of  oscillation  over  toe  range. 


11 


It  can  be  shown  that  the  approximation  holds,  provided  the  range  interval  is  large  enough  to  admit  many 
oscillations,  or  small  enough  so  that  the  transfer  function  does  not  change  appreciably.  If  the  range 
oscillations  behave  such  that  the  wavelength  is  equal  to  the  familiar  quotient  of  sound  speed  and 
frequency,  the  range  of  integration  should  extend  over  the  1/3-octave  range  intervals.  The  rms  average 
over  the  1/3-octave  band  can  be  replaced  by  an  average  over  the  1/3-octave  range  interval.  Note  that  the 
wavelength  of  the  oscillations  in  range  may  not  obey  the  relation  above,  but  it  is  not  a  strict  requirement 
for  the  argument  to  hold. 

The  predictions  from  the  FFP  were  range-averaged  in  an  attempt  to  remove  the  unwanted 
oscillations  in  the  Lookup  data  base  files.  The  arithmetic  standard  deviations  of  the  levels  over  these 
range  intervals  were  also  calculated  and  stored.  The  standard  deviation  is  a  useful  indicator  of  the  possible 
variation  in  levels  due  to  small  changes  in  position,  frequency,  or  environmental  conditions.  Future 
improvements  to  the  Lookup  program  might  involve  the  use  of  this  quantity  for  a  prediction  of  the 
anticipated  variation  in  levels,  or  the  probability  of  exceeding  certain  thresholds. 

Occasionally  the  FFP  yielded  predictions  at  so  few  ranges  that  some  range  intervals  contained  no 
prediction.  In  cases  where  there  are  no  predicted  levels,  the  data  base  entries  were  set  equal  to  100.0  dB 
and  the  corresponding  standard  deviations  were  set  equal  to  -1.0  dB. 

In  the  initial  release  of  the  data  base,  several  atmospheric  and  environmental  parameters  varied  in 
the  table.  Every  combination  of  the  following  parameters  was  used:  hard  and  porous  ground:  linear  and 
logarithmic  sound  speed  profiles;  profile  speeds  ranging  from  about  -12  m/s  to  about  12  m/s  in  steps  of. 
2  m/s;  source  heights  of  2  m,  5  m,  15  m,  and  100  m;  and  frequencies  from  1  Hz  to  1600  Hz,  on 
1/3-octave  centers.  In  all,  a  total  of  370  separate  cases  were  prepared  for  the  table. 

Though  this  table  does  not  (cannot)  cover  all  possible  situations  that  might  occur  in  practice,  it  does 
provide  a  starting  place  for  near-ground  propagation  to  medium  range  (up  to  about  3  km).  Beyond  3  km 
or  so,  other  influences  often  become  more  important,  such  as  hilly  terrain  or  turbulence.  The  nature  of 
the  atmosphere  well  above  conventional  surface-based  meteorological  sensors  has  a  more  pronounced 
effect  on  levels  at  longer  ranges.  Nevertheless,  computations  for  these  cases  were  performed,  and  stored 
in  the  table  for  distances  up  to  20  km  in  range. 
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SUMMARY 


This  study  has  created  a  method  to  quickly  predict  sound  propagation  from  outside  sources  through  the 
Lookup  program  and  data  base.  Based  on  the  assumption  that  many  situations  in  outdoor  propagation  tend  to 
recur,  the  Lookup  program  speeds  the  prediction  process  by  accessing  a  data  base  of  information  describing 
typical  situations,  and  using  that  information  to  estimate  received  sound  levels. 

The  user  provides  the  Lookup  program  with  an  input  file  that  contains  the  physieal  parameters  of  the 
propagation  problem,  e.g.,  the  receiver  heights,  winds  speed  and  direction,  type  of  ground  surface,  and  so  forth. 
The  Lookup  program  estimates  sound  propagation  by  matching  the  entered  information  against  its  data  base  of 
typical  parameters,  and  doing  a  calculation  based  on  the  closest  match. 
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APPENDIX  A:  Lookup  Program  Source  Code 


Use  lookup  table  to  estimate  sound  propagation. 


CHARACTER  Dir*80 , Best*80, Index*80 
INTEGER  Iband, Prof , Rexs , Rold, Rint 
LOGICAL  Fend,Fopen 

REAL  Srchi , Rechi , Range, Azimth, Grclas , Ruf hi 
REAL  Tchi (100), Tc (100), Rhhi (100),Rh(100) 

REAL  Whi (100) ,Wx(100) ,Wy(100) 

REAL  TobSrc (0:43), Ref dis , Usrcrv ( 0 : 43 ) , TobFlt (0:43) 

REAL  Tc  0  ,  Dt dh ,  RhO  ,  S 1 0  ,  W1 0  ,  W1  Opro ,  Fr ,  Airab 

REAL  Delair ( 0 : 43 ) , Delstp ( 0 : 43 ) , Delexs ( 0 : 43 ) , Delold ( 0 : 43 ) 

REAL  R, Rng, L, A, C, U, Delspr , Delint ( 0 : 43 ) , R2 , S, XI , X2 
REAL  TobAwt (0:43) , TobCwt (0:43) , TobUwt (0:43) , Ovral 1 

EXTERNAL  Airab, Ovrall 

- Do  all  the  parameter  reading. 

CALL  Rdver ( ) 

CALL  Rdname (Dir, Index) 

CALL  Rdgeom ( Srchi , Rechi , Range , Azimth) 

CALL  Rdgrnd (Grclas, Rufhi) 

CALL  Rdtc(Tchi,Tc) 

CALL  Rdrh(Rhhi,Rh) 

CALL  Rdwswd(Whi,Wx,Wy) 

CALL  Rdlvls (TobSrc, Refdis, Usrcrv) 

- Fit  met  profiles;  choose  speed  profile;  select  "best"  file. 

CALL  Fittc(Tchi,Tc,TcO,Dtdh) 

CALL  Fitrh (Rhhi, Rh, RhO) 

CALL  Fitw(Whi,Wx,Wy, Rufhi, Azimth, W10,W10pro) 

CALL  Chprof (Dtdh,W10,W10pro, Prof ,S10) 

CALL  Fndfil (Prof , SIO, Grclas, Srchi, Rechi, Dir, Index, Best) 

- Molecular  absorption  spectrum  (dB/m)  offset  from  STP . 

DO  10  Iband=0,43 

Fr  =  10.0**(Iband/10.0) 

Delair ( Iband)  =  8 . 686*Airab (Fr , TcO , RhO ,  1 . 0 ) 

Delstp ( Iband)  =  8 . 686*Airab (Fr , 15 . 0 , 70 . 0 , 1 . 0 ) 

10  CONTINUE 

- Prepare  to  read  the  excess  levels  file. 

Fopen  =  .FALSE. 

Fend  =  .FALSE. 

Rexs  =  -1 
Rold  =  -1 

WRITE  (*,*)  'BEGIN  RANGE,  FLAT,  A,  C  AND  USER' 

DO  70  Rint=0,43 

- Choose  the  range  number  for  interpolation  or  sweep. 

IF  (Range .GT. 0 . 0)  THEN 
R  =  10 . 0*ALOG10 (Range) 

ELSE 

R  =  FLOAT (Rint) 

ENDIF 

- Fetch  the  next  good  spectrum. 

30  IF  ( .NOT. Fend. AND. Rexs .LT.R)  THEN 

- Save  the  old  spectrum  first. 

Rold  =  Rexs 
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c 

c 


c 

c 


c 

c 


c 

c 


c 

c 

c 

c 


c 

c 

c 

c 


c 

c 


c 

40 


CALL  Hold (Delexs, Delold) 

- Here  we  attempt  to  get  the  new  one  (may  be  EOF,  etc.) . 

CALL  Rdexs (Rexs , Delexs , Fopen, Fend, Dir , Best ) 

GOTO  3  0 

- Save  the  good  one  the  last  trip  through  (reached  EOF) . 

ELSE  IF  (Fend)  THEN 

CALL  Hold(Delold, Delexs) 

ENDIF 


- See  if  interpolation  is  possible. 

IF  ( (Range. GT. 0.0) .AND. (Rold.GE. 0.0) .AND. ( .NOT. Fend) )  THEN 

CALL  NtrplS (FLOAT(Rold) , Delold, FLOAT (Rexs ) , Delexs , R, Delint ) 
CALL  Hold (Delint, Delexs) 

ENDIF 


- Spherical  distance. 

Rng  =  10.0**(R/10.0) 

R2  =  SQRT( (Srchi-Rechi) **2+Rng**2) 

Delspr  =  20 . 0*ALOG10 (R2/Refdis) 

- Apply  spherical  spreading  and  molec.  absorption. 

DO  40  Iband=0,43 

- Source  level. 

S  =  TobSrc (Iband) 

TobFlt (Iband)  =  S 

- Spherical  spreading  plus  air  absorb. 

XI  =  -Delspr-Delair (Iband) *R2 

- Remove  std.  air  from  table  attn.  Cut  losses  to  -25. 

X2  =  AMAXl (Delexs (Iband) +Delair (Iband) *R2, -25 . 0) 

- Only  modify  levels  if  source  band  is  set.  Stop  at  -100. 

■  IF  (S .GT. -100 . 0 )  TobFlt ( Iband)  =  AMAXl (S+X1+X2, -100 . 0) 

CONTINUE 


C  - Apply  flat  (none).  A,  C,  and  User  weighting. 

C ALL  F rwt ng ( TobF 1 1 , TobAwt , TobCwt , TobUwt , Us  r c  r v ) 

L  =  Ovrall (TobFlt) 

A  =  Ovrall (TobAwt ) 

C  =  Ovrall (TobCwt) 

U  =  Ovrall (TobUwt) 

C 

0  - Output  the  Range,  Sum,  A,  C,  and  User  SEL  dB . 

WRITE  (*,50)  Rng,L,A,C,U 
50  FORMATdX,  F7 . 1, 4  (IX,  F6.1)  ) 

C 

C  - Print  the  interpolated  spectrum. 

IF  (Range. GT. 0.0)  THEN 

WRITE  (*,*)  'RECEIVER  SPECTRUM,  1/3-OCTAVE  BANDS  0-43 
WRITE  (*,60)  (Tobf It (Iband) , lband=0, 43) 

60  FORMAT  (5F7.1) 

GOTO  80 
ENDIF 
70  CONTINUE 
80  END 

REAL  FUNCTION  Airab (Fr , Tc, Rh, P) 

C - - - 

C  Absorption  of  sound  in  air  by  molecular  processes. 

C  [1]  American  National  Standard  Method  for  the  Calculation  of  the 

C  Absorption  of  Sound  by  the  Atmosphere,  ANSI  SI. 26-1978. 
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C  [2]  H.  E.  Bass,  L.  C.  Sutherland,  J.  Piercy  and  L.  Evans, 

C  "Absorption  of  Sound  by  the  Atmosphere,"  in  'Physical 

C  Acoustics'  (1984). 

C - - - 

C  Fr  --  Frequency  (Hz) 

C  Tc  --  Temperature  (deg  C) 

C  Rh  --  Relative  Humidity  (%)  ( 0  <  Rh  <  100) 

C  p  __  Pressure  (standard  sea-level  atm) 

C - - - 

REAL  Alpha, Fr, Frn, Fro, H, P, Psatpo, Rh, T, TO , T20ot , Tc 

- Temperature  in  Kelvins,  and  reference  temperature  20  deg  C. 

T  =  Tc+273.15 
TO  =  293 . 15 
T20ot  =  TO/T 

- Saturation  pressure  ([2],  Eqn.  72). 

Psatpo  =  10. 0**(8. 422-10. 05916*T20ot+5. 023*AL0G10(T20ot) ) 

- Percent  mole  fraction  of  water  vapor  ([1],  Eqn.  DIO). 

H  =  Rh*Psatpo/P 

- Oxygen  and  nitrogen  relaxation  frequencies  ([1],  Eqs .  8  &  9). 

Fro  =  P*(24.0+4.41E4*H*(0.05+H)/(0.391+H) ) 

Frn  =  P*SQRT(T20ot) * (9 . 0  +  3 50 . 0*H*EXP ( -6 . 142  * ( T20ot**0 . 3 33 3 -1 . 0) ) ) 

- Absorption  coefficient.  Alpha  (Nepers/m)  ([1],  Eqn.  10). 

Alpha  =  0.01278*EXP(-2239.1/T)/(Fro+Fr*Fr/Fro) 

Alpha  =  0.10690*EXP(-3352.0/T)/(Frn+Fr*Fr/Frn)  +  Alpha 
Alpha  =  Fr*Fr*(1.84E-ll/SQRT(T20ot)/P+T20ot**2.5*Alpha) 

Airab  =  Alpha 

RETURN 
END 

SUBROUTINE  Chprof (Dtdh, W10,W10pro, Prof ,  SIO) 


Choose  the  sound  speed  profile  type  from  met  profiles 

-  linear:  Prof=0 

-  log:  Prof=l 


INTEGER  Prof 

REAL  Dtdh, SIO, W10,W10pro 

- Decide  which  profile  type  to  use. 

IF  (WIO.GT.2.0)  THEN 

- If  wind  magnitude  >  2  m/s  use  log  speed  profile. 

Prof  =  1 
SIO  =  WlOpro 

ELSE 

- Use  a  linear  profile,  based  on  the  temperature. 

Prof  =  0 

SIO  =  Dtdh*0. 61*10.0 
END  IF 

RETURN 

END 

REAL  FUNCTION  Dbsum(X,Y} 


Add  energies  dB (E (X) +E (Y) )  of  two  uncorrelated  signals. 


REAL  Lo,Hi,X,Y 

- Lower  and  higher  values  of  X  and  Y  to  avoid  overflow. 
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Lo  =  AMIN1(X,Y) 

Hi  =  AMAXl (X, Y) 

IF -(Lo.LE. -100.0)  THEN  ^ 

- Don't  add  in  values  of  0.0,  (zero  energy  flag)  . 

Dbsum  =  AMAXl ( Hi , -10 0 . 0 ) 

ELSE  ,  1  4-  jT, 

r  - Convert  to  energy,  add,  and  convert  back  to  dB . 

C  Dbsum  =  10.0*ALOGl0(10.0**(Hi/10.0)+10.0**(Lo/10.0) ) 

Dbsum  =  Hi+10.0*ALOGl0{1.0+10.0**((Lo-Hi)/10.0)) 

END  IF 
C 

RETURN 
F  D 

INTEGER  FUNCTION  Endstr(C) 

C  Function  returns  the  position  of  the  last  non-blank^charater^^^ 

CHARACTER  C*80 
INTEGER  Last  ■ 

C 

DO  10  Last=80,l,-1 

IF  (C (Last : Last) .NE. '  ')  GOTO  20 

10  CONTINUE 
C 

20  Endstr  =  Last 
RETURN 
END 

SUBROUTINE  Fitrh (Rhhi , Rh, RhO ) 

C  Find  the  relative  humidity  nearest  the  surface  from^the ^profile 

INTEGER  I,N 

REAL  Low, Rhhi ( 100 ), Rh ( 100) , RhO 

- Count  up  the  number  of  data  points. 

DO  10  N=0,99 

IF  (Rhhi (N+1) .LT. 0 . 0)  GOTO  20 
10  CONTINUE 

20  IF  (N.EQ. 0)  THEN 

2  - Assume  standard  humidity . 

RhO  =70.0 

ELSE  IF  (N.EQ.l)  THEN 

2  - Use  the  only  available  value. 

RhO  =  Rh(l) 

ELSE 

2  - Use  value  closest  to  the  surface. 

Low  =  Rhhi ( 1 ) 

RhO  =  Rh(l) 

DO  30  1=2, N 

IF  (Rhhi ( I) .LT.Low)  THEN 
Low  =  Rhhi ( I ) 

RhO  =  Rh(I) 

ENDIF 

3  0  CONTINUE 

ENDIF 
C 

RETURN 

END 

SUBROUTINE  Fittc (Tchi , Tc , TcO , Dtdh) 

C  Find  the  slope  and  intercept  of  the  temperature  profile^ _ 

REAL  Dtdh,  TcO,  Tc  (100)  ,  Tchi  (100) 

INTEGER  N 


C 

C 
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- Count  up  the  number  of  data  points. 

DO  10  N=0,99 

IF  (Tchi(N+l) .LT.0.0)  GOTO  20 
10  CONTINUE 

20  IF  (N.EQ.O)  THEN 

- Assume  standard  temp,  lapse  rate. 

Dtdh  =  -0.0098 
TcO  =15.0 

ELSE  IF  (N.EQ.l)  THEN 

- Assume  lapse  rate,  find  To  at  the  ground. 

Dtdh  =  -0.0098 
TcO  =  Tc (1) -Tchi (1) *Dtdh 
ELSE 

- Use  slope  and  intercept  of  best-fit  line. 

CALL  Linf it (Tchi , Tc, N, TcO , Dtdh) 

ENDIF 

RETURN 

END 

SUBROUTINE  Fitw (Whi , Wx, Wy, Ruf hi , Azimth, WIO , WlOpro) 


Find  the  best  log  wind  profile. 


REAL  Azimth, Pi , Rufhi , Yx, Yy 

REAL  W10,W10pro,W10x,W10y,Whi (100) ,Whilg(100) ,Wx(100) ,Wy(100) 

.  INTEGER  I,N 

- Figure  out  Pi . 

Pi  =  4.0*ATAN(1.0) 

- Count  up  the  number  of  data  points. 

DO  10  N=0,99 

IF  (Whi (N+1) .LT. 0 . 0)  GOTO  20 
10  CONTINUE 

20  IF  (N.EQ.O)  THEN 

- Default  wind  speed  at  2 . 0  m/s  along  propagation  azimuth. 

WIO  =  2.0 
WlOpro  =2.0 
ELSE 

- Make  a  loglO  height  array. 

DO  30  1=1, N 

Whilg(I)  =  ALOGIO (Whi (I) ) 

3  0  CONTINUE 

- Insert  false  datum  at  the  roughness  height. 

Whi Ig (N+1)  =  ALOGIO (Rufhi) 

Wx(N+l)  =0.0 
Wy(N+l)  =0.0 

- Least  squares  fit  on  (Whilg [ ] , Wx [ ] )  and  (Whilg [ ] , Wy [ ] ) . 

CALL  Linfit (Whilg, Wx, N+1 , Yx, WlOx) 

CALL  Linfit (Whilg, Wy, N+1 , Yy, WlOy ) 

- Wind  magnitude. 

WIO  =  SQRT(W10x**2+W10y**2) 

-  Wind  flow  along  the  propagation  azimuth. 

-  WlOpro  =  (WlOx  +  WlOy)  (dot)  (sin (Azimth) x  +  cos (Azimth)y) 

WlOpro  =  W10x*SIN(Azimth*Pi/180 . 0) +W10y*COS (Azimth*Pi/180 . 0) 
ENDIF 

RETURN 
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END 

SUBROUTINE  Fndf il (Tprf , Tspd, Tgrd, Tsrc, Tree, Dir, Index, Bf il) 


c--- 

c 

Find  the  best  file  for  this  propagation  condition. 

c — 

c 

B  -  Best  match. 

c 

N  -  Newest  match. 

c 

c--- 

T  -  Target  match. 

CHARACTER  Dir*80, Index*80 , Bf il*80 , Nf il*80 
INTEGER  Ends t r, Nprf, Tprf 
LOGICAL  Match 

REAL  Bspd, Nspd, Tspd, Difspd,  Bgrd, Ngrd, Tgrd, Difgrd 
REAL  Bsrc, Nsrc, Tsrc, Dif sre,  Brec , Nrec , Tree , Dif rec ,  Tol 
EXTERNAL  Endstr 
C 

Match  =  .FALSE. 

Tol  =  0.001 
C 

OPEN  ( 10 , FILE=Dir ( 1 :Endstr (Dir) ) / /Index, STATUS =' OLD ' ,ERR=30) 
10  READ  (10, *,ERR=30,END=20)  Nf i 1 , Nprf , Nspd, Ngrd, Nsrc , Nrec 

- Skip  if  not  target  profile  type. 

IF  (Nprf .EQ. Tprf )  THEN 

- If  not  yet  a  match,  set  these  best  parameters. 

IF  (.NOT. Match)  THEN 
Bspd  =  Nspd 
Bgrd  =  Ngrd 
Bsrc  =  Nsrc 
Brec  =  Nrec 
Bfil  =  Nfil 
Match  =  .TRUE. 

ELSE 

- There  is  at  least  one  match;  fit  best  speed. 

Difspd  =  ABS (Nspd-Tspd) -ABS (Bspd-Tspd) 

IF  (Difspd.LT.0.0)  THEN 
Bspd  =  Nspd 
Bgrd  =  Ngrd 
Bsrc  =  Nsrc 
Brec  =  Nrec 
Bfil  =  Nfil 

ELSE  IF  (Difspd.LT.Tol)  THEN 

- Same  speed;  fit  best  ground  type. 

Difgrd  =  ABS (Ngrd-Tgrd) -ABS (Bgrd-Tgrd) 

IF  (Difgrd. LT. 0 . 0)  THEN 
Bgrd  =  Ngrd 
Bsrc  =  Nsrc 
Brec  =  Nrec 
Bfil  =  Nfil 

ELSE  IF  (Difgrd. LT. Tol)  THEN 

- Same  speed, ground;  fit  best  source  height. 

Difsrc  =  ABS (Nsrc-Tsrc) -ABS (Bsrc-Tsrc) 

IF  (Difsrc.LT.0.0)  THEN 
Bsrc  =  Nsrc 
Brec  =  Nrec 
Bfil  =  Nfil 

ELSE  IF  (Difsrc.LT.Tol)  THEN 

- Same  speed,  ground,  source;,  fit  best  rec. 

Difrec  =  ABS (Nrec-Trec) -ABS (Brec-Trec) 

IF  (Difrec .LT. 0 . 0 )  THEN 
Brec  -  Nrec 
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Bfil  =  Nfil 
ENDIF 
END  IF 
ENDIF 
ENDIF 
ENDIF 
ENDIF 
GOTO  10 


- Flag  a  missing  decay  table  as  a  fatal  error. 

0  IF  (.NOT. Match)  THEN 

PRINT  *, 'ERROR:  No  decay  table  match  found.' 

WRITE  (*,*)  ' INDEX  FILE  NAME: ', Dir (1 :Endstr (Dir) ) //Index 

STOP 

ENDIF 

RETURN 


-Some 
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(* 

★ 
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★ 
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★ 
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WRITE 
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★ 

WRITE 

(* 

★ 

WRITE 

(* 

★ 

WRITE 

(* 

★ 

WRITE 

(* 

* 

STOP 

END 

kind 

) 

) 

) 


of  reading  error. 

'ERROR:  bad  entry  in  index  file.' 

'INDEX  FILE  NAME: ', Dir (l:Endstr( Dir) ) //Index 
'expecting  file  name  with  parameters,' 
'instead  found:' 

'FILE  NAME: ' ,Nfil 
'PROFILE  TyPE:',Nprf 
'PROFILE  SPEED: ',Nspd 
'  GROUND  CLASS : ' , Ngrd 
'SOURCE  HEIGHT: ' ,Nsrc 
'RECEIVER  HEIGHT: ' , Nrec 


SUBROUTINE  Frwtng (Tobin, TobAwt , TobCwt , TobUwt , Usrcrv) 

Apply  frequency  weighting  curve  to  1/3-octave  bands  0  to  43 


INTEGER  Iband 

REAL  Acrv(0:43) ,Ccrv(0:43) ,Usrcrv(0:43) 

REAL  Tobin (0:43), TobAwt (0:43), TobCwt (0:43), TobUwt (0:43) 


DATA  Acrv/ 


> 

-148.5, 

-140.5, 

> 

-108.7, 

-100.8, 

> 

-70.4, 

-63.3, 

> 

-39.4, 

-34.6, 

> 

-19.1, 

-16.0, 

> 

-6.6, 

-4.8, 

> 

0.0, 

0.5, 

> 

1.1, 

0.9, 

> 

-2.4, 

-4.3, 

132.6, 

-124.6,  - 

116.6, 

-93.0, 

-85.3, 

-77.7, 

-56.6, 

-50.4, 

-44.7, 

-30.2, 

-26.1, 

-22.5, 

-13.3, 

-10.8, 

-8.6, 

-3.2, 

-1.9, 

-0.8, 

0.9, 

1.1, 

1.2, 

0.5, 

-0.1, 

-1.1, 

-6.6, 

-9.3/ 

10 


DATA  Ccrv/ 


> 

-52.5, 

-48.5, 

-44.5, 

-40.5, 

-36.6, 

> 

-32.6, 

-28.8, 

-24.9, 

-21.2, 

-17.6, 

> 

-14.3, 

-11.2, 

-8.5, 

-6.2, 

-4.4, 

> 

-3.0, 

-1.9, 

-1.2, 

-0.8, 

-0.5, 

> 

-0.2, 

-0.1, 

0.0, 

0.0, 

0.0, 

> 

0.0, 

0.0, 

0.0, 

0.0, 

0.0, 

> 

0.0, 

0.0, 

0.0, 

-0.1, 

-0.2, 

> 

-0.5, 

-0.8, 

-1.2, 

-1.9, 

-3.0, 

> 

-4.4, 

-6.2, 

-8.5, 

-11.2/ 

- Apply  weightings  to  the  spectra. 

DO  10  Iband=0,43 

TobAwt (Iband)  =  AMAXl ( -10 0 . 0 , Tobin ( Iband) + Acrv ( Iband) ) 
TobCwt (Iband)  =  AMAXl ( -100 . 0 , Tobin ( Iband) +Ccrv ( Iband) ) 
TobUwt(Iband)  =  AMAXl  ( -100 . 0 , Tobin ( Iband) -i-Usrcrv ( Iband)  ) 
CONTINUE 
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c 

RETURN 

END 

SUBROUTINE  Hold(X,Y) 

C- - 

C  Copy  one  spectrum  X  into  spectrum  Y. 

C - 

INTEGER  Iband 
REAL  X(0:43) ,Y(0:43) 

C 

DO  10  lband=0 , 43 

Y( Iband)  =  X( Iband) 

10  CONTINUE 
C 

RETURN 

END 

SUBROUTINE  Linf it (X, Y, Ndata, YO , Slope) 


C  Least  squares  linear  fit  y=Slope*x+Y0  to  data  y ( i ) , x ( i ) , i-1 , Ndata^ 

C  [1]  Press,  et  al . ,  Numerical . Recipes ,  Cambridge,  pp .  504-509. 

INTEGER  I , Ndata 

REAL  X(Ndata) ,Y(Ndata) ,Y0, Slope,  Sxon , Sx, Sy , St2 , T 
C 

Sx  =  0 . 0 
Sy  =  0.0 
St2  =0.0 
Slope  =0.0 
C 

Q  - Sum  over  X  and  over  Y . 

DO  10  1=1, Ndata 
Sx  =  Sx+X(I) 

Sy  =  Sy+Y(I) 

10  CONTINUE 

Sxon  =  Sx/FLOAT (Ndata) 

C 

C  - Sum  of  squared  differences. 

DO  20  1=1, Ndata 
T  =  X(I) -Sxon 
St2  =  St2+T*T 
Slope  =  Slope+T*Y(I) 

20  CONTINUE 
C 

C  - Intercept  and  slope. 

Slope  =  Slope/St2 

YO  =  (Sy-Sx*Slope) /FLOAT (Ndata) 

C 

RETURN 

END 

SUBROUTINE  NtrplS ( A, X, B , Y ,  C ,  Z ) 

- - - 

C  Interpolates  between  spectra  X  and  Y,  based  on  value  of  C  in  (A,B) 

C - 

INTEGER  Iband 

REAL  Ratio,A,B,C,X(0:43) ,Y(0:43)  ,Z(0:43) 

C 

C  - Find  the  fraction  of  C  of  the  way  between  A  and  B. 

Ratio  =  (C-A) / (B-A) 

C 

C  - - Let's  just  don't  extrapolate. 

IF  (Ratio .GT. 1 . 0 .OR.Ratio.LT. 0 . 0)  THEN 

WRITE  (*,*)  'ERROR:  interpolation  attempted  outside  interval' 
WRITE  (*,*)  'looking  for:',C,'in  interval :', CMPLX (A, B) 

STOP 
END  IF 
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c 

c 


10 

c 


c — 
c 

c - 


c 

c 


10 

c 

c 


c - 

c 

c - 


c 

c 


10 

c 

c 


c 

c 

20 

30 

c 

c 

40 

50 

C 


60 

C 

C 

C 

C 

C 

C 


- Do  all  of  the  bands. 

DO  10  Iband=0,43 

Z(Iband)  =X(Iband)  +  (Y ( Iband) -X ( Iband) ) *Ratio 
CONTINUE 

RETURN 

END 

REAL  FUNCTION  Ovrall(Tob) 


Sum  over  the  1/3 -octave  bands  0  to  43. 

INTEGER  Iband 
REAL  Tob(0:43),  Leq 
REAL  Dbsum 
EXTERNAL  Dbsum 

- Find  the  overall  level. 

Leq  =  -100 . 0 
DO  10  Iband=l,43 

Leq  =  Dbsum (Tob( Iband) ,Leq) 
CONTINUE 


Ovrall  =  Leq 

RETURN 

END 

SUBROUTINE  Patch (X, Bad) 

Fill  outside  band  edge  and  interpolate  to  fill  holes^in^spectrum^ 

INTEGER  Iband, Jband, First, Last 

LOGICAL  Bad 

REAL  X( 0 : 43 ), Ratio 

Bad  =  . FALSE . 

- Find  the  first  good  band. 

DO  10  First=0,43 

IF  (X(First) .LT.lOO .0)  GOTO  20 
CONTINUE 

- All  bands  are  bad. 

Bad  :=  .TRUE. 

RETURN 

- Find  the  last  good  band. 

DO  30  Last=43,0,-1 

IF  (X(Last) .LT. 100 . 0)  GOTO  40 
CONTINUE 

_ Fill  all  bands  on  the  outer  edges. 

DO  50  Iband=First-l, 0, -1 
X(Iband)  =  X(First) 

CONTINUE 

DO  60  Iband=Last+l,  43 
X(Iband)  =  X(Last) 

CONTINUE 

- Scan  for  holes  in  the  spectrum. 

DO  90  Iband=First+l, Last-1 

- Check  if  missing  point. 

IF  (X(Iband) .GE.lOO .0)  THEN 

- Find  the  next  good  band  (last  band  is  OK) . 
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70 

C 

C 

80 


C 

90 

C 


C - 

C 

C - 


C 

C 


c 

c 

1 

5 

C 

C 


C 

C 


C 

C 

10 


c 

c 

20 


C 

C 

40 


C-- 

C 


DO  70  Jband=Iband+l , Last 

IF  (X(Jband) .LT.lOO .0)  GOTO  80 
CONTINUE 

_ Use  the  Jband  and  Iband-1  to  interpolate  on. 

Ratio  =  1 . 0/FLOAT( Jband- Iband+1) 

X(Iband)  =  X(Iband-l) + (X(Jband) -X(Iband-l) ) *Ratio 

ENDIF 


CONTINUE 


RETURN 

SUBROUTINE  Rdexs (Rnum, Exsdb, Fopen, Fend, Dir , Best) 


Read  the  excess  levels  file. 

CHARACTER  Best*80 , Dir*80 , Header*20 
INTEGER  Iband, Endstr , Rnum 
LOGICAL  Fend , Fopen , Bad 
REAL  Exsdb (0:43) 

EXTERNAL  Endstr 


- __Take  care  of  opening  the  "Best"  file. 

IF  (.NOT. Fopen)  THEN 

OPEN  (20 , File=Dir ( 1 :Endstr (Dir) ) / /Best , STATUS= 
Fend  =  .FALSE. 

Fopen  =  . TRUE . 

REWIND  (UNIT=20) 

ENDIF 


'OLD' ,ERR=10) 


- Read  the  range . 

READ  (20, 5,ERR=20,END=40)  Header 
FORMAT  (A20) 

READ  (20,*)  Rnum 


- Read  average  spectrum 

READ  (20, *,END=20) 

READ  (20, *,ERR=20,END=20) 


,  bands  0-4,  5-9,  etc.  through 
(Exsdb ( Iband) , lband=0 , 43 ) 


43  . 


_ Try  and  patch  empty  holes  in  the  decay  spectra. 

CALL  Patch (Exsdb, Bad) 

IF  (Bad)  GOTO  1 
RETURN 

- Bad  file  name  or  missing  file  error. 

WRITE  (*,*)  'ERROR;  decay  file  doesn''t  open:' 

WRITE  (*,*)  ' Directory/File :', Dir ( 1 : Endstr (Dir ))/ /Best 

STOP 

- VJe  found  end  of  file  too  early,  or  bad  numbers.  ^ 

WRITE  (*,*)  'ERROR:  early  end  of  file  or  bad  spectrum: 

WRITE  (*,*)  'Directory/File :', Dir ( 1 :Endstr (Dir) ) //Best 

WRITE  (  * , * )  ' Range  number : ' , Rnum 

WRITE  (*,*)  'Spectrum:' 

WRITE  (*,30)  (Exsdb ( Iband) , lband=0 , 43 ) 

FORMAT  (5F7.1) 

STOP 

- Found  end  of  file. 

Fend  =  .TRUE. 

RETURN 

END  .  . 

SUBROUTINE  Rdgeom(Srchi , Rechi , Range , Azimth) 


Read  in  the  geometry  of  the  propagation  problem. 
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- - 

CHARACTER  Header *20 
REAL  Azimth, Range, Srchi, Rechi 
10  FORMAT  (A20) 

- Source  height  (m) . 

READ  (*, 10,ERR=30,END=20)  Header 
IF  (Header (1: 13 ) .NE. 'SOURCE  HEIGHT')  GOTO  30 
READ  (*, *,ERR=30,END=20)  Srchi 

- Reciever  height  (m) . 

READ  (*, 10,ERR=30,END=20)  Header 
IF  (Header (1:15) .NE. 'RECEIVER  HEIGHT')  GOTO  30 
READ  (*, *,ERR=30,END=20)  Rechi 
IF  (Srchi .LT. 0 . 0 .OR. Rechi .LT. 0 . 0)  THEN 

WRITE  (*,*)  'ERROR:  Source  and  receiver  must  be  above  ground.' 
WRITE  (*,*)  'Instead,  f ound :', Srchi , Rechi 
STOP 
END  IF 

- Range  for  the  prediction  (m)  (0=ALL). 

READ  (*, 10,ERR=30,END=20)  Header 
IF  (Header (1: 14) .NE. 'RECEIVER  RANGE' )  GOTO  30 
READ  (*, *,ERR=30,END=20)  Range 

- Azimuth  for  propagation  (deg) . 

READ  (*, 10,ERR=30,END=20)  Header 

IF  (Header (1: 16) .NE. 'RECEIVER  AZIMUTH')  GOTO  30 
READ  (*, *,ERR=30,END=20)  Azimth 
RETURN 

20  Header  =  'EOF' 

30  WRITE  (*,40)  Srchi, Rechi, Range, Azimth,  Header 

40  FORMAT (' ERROR :  expecting  a  header,  followed  by  a  value:'/ 

>  5X, 'SOURCE  HEIGHT'/  F15.2/  5X, 'RECEIVER  HEIGHT'/  F15.2/ 

>  5X, 'RECEIVER  RANGE ' /  F15 . 2 /  5X, 'RECEIVER  AZIMUTH'/  F15.2/ 

>  'error  found  with:',A20) 

STOP 

END 

SUBROUTINE  Rdgrnd (Grclas , Rufhi ) 


Read  in  the  ground  conditions. 


CHARACTER  Header *20 
REAL  Grclas, Rufhi 
10  FORMAT  (A20) 

- Ground  classification  type  (hard=0,  soft=l). 

READ  (*, 10,ERR=30,END=20)  Header 
IF  (Header (1: 12) .NE. 'GROUND  CLASS')  GOTO  30 
READ  (*, *,ERR=30,END=20)  Grclas 
IF  (Grclas. LT. 0.0. OR. Grclas. GT. 1.0)  THEN 

WRITE  (*,*)  'ERROR:  Ground  class  must  be  between  0.0  and  1.0.' 
WRITE  (*,*)  'Instead,  found :', Grclas 
STOP 
ENDIF 

- Ground  roughness  height  (m)  . 

READ  (*, 10,ERR=30,END=20)  Header 
IF  (Header (1: 12) .NE. 'GROUND  ROUGH')  GOTO  30 
READ  (*, *,ERR=30,END=20)  Rufhi 
IF  (Rufhi.LE.0.0)  THEN 

WRITE  (*,*)  'ERROR:  Ground  roughness  height  must  exceed  0.0.' 
WRITE  (*,*)  'Instead  f ound :', Rufhi 
STOP 
ENDIF 


RETURN 


20 

30 
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C 

C 

10 


C 

c 


c 

c 


c 

c 

20 

30 


40 


C 

C 

50 

60 


C 

C 

70 

80 


C 


C-- 

C 

C-- 

C 

C 


Header  =  'EOF' 

WRITE  (*,40)  Grclas,Rufhi, Header 

FORMAT (' ERROR :  expecting  a  header,  followed  by  a  value:'/ 

>  5X, 'GROUND  CLASS'/  F15.2/  5X, 'GROUND  ROUGH'/  F15.2/ 

>  'error  found  with:',A20) 

STOP 

END 

SUBROUTINE  Rdlvls ( TobSrc , Ref dis , Usrcrv) 


Calls  routines  to  read  source  spectrum  and  weighting. 


CHARACTER  Header*20 
INTEGER  Iband 

REAL  Refdis, TobSrc (0:43) , Usrcrv ( 0 : 43 ) 


- -Source  spectrum. 

READ  (*, 10,ERR=30,END=20)  Header 
FORMAT  (A20) 

IF  (Header (1:18) .NE. 'SOURCE  LEVEL  SPECT' )  GOTO  30 
READ  (*, *,ERR=30,END=20)  (TobSrc ( Iband) , lband=0 , 43 ) 


- Reference  distance. 

READ  (*, 10,ERR=60,END=50)  Header 
IF  (Header (1: 14) .NE. 'REFERENCE  DIST')  GOTO  60 
READ  (*, *,ERR=60,END=50)  Refdis 
IF  (Refdis.LE.0.0)  GOTO  60 


- -User  weighting  spectrum. 

READ  (*, 10,ERR=80,END=70)  Header 

IF  (Header (1:20) .NE. 'USER  WEIGHTING  SPECT')  GOTO  80 
READ  (*, *,ERR=80,END=70)  (Usrcrv ( Iband) , lband=0 , 43 ) 
RETURN 


- Error  in  source  spectrum. 

Header  =  'EOF' 

WRITE  (*,*)  'ERROR:  expecting  source  spectrum.' 
WRITE  (*,*)  'Instead,  found: ', Header 
WRITE  (*,*)  'Spectrum  =' 

WRITE  (*,40)  (TobSrc ( Iband) , lband=0 , 43 ) 

FORMAT  (5F7.1) 

STOP 


- Error  in  reference  distance. 

Header  =  'EOF'  ^ 

WRITE  (*,*)  'ERROR:  expecting  positive  reference  distance. 
WRITE  (*,*)  'Instead,  f ound :', Header , '  Ref dis= ' , Ref dis 
STOP 


- Error  in  user  weighting  spectrum. 

Header  =  'EOF'  ^  ^  ^ 

WRITE  (*,*)  'ERROR:  expecting  user  weighting  spectrum. 

WRITE  (*,*)  'Instead,  f ound :', Header 
WRITE  (*,*)  'Spectrum  := ' 

WRITE  (*,40)  (Usrcrv ( Iband) , lband=0 , 43 ) 

STOP 

END 

SUBROUTINE  Rdname ( Dir , Index) 

Reads  the  directory  name  and  index  file  name  from  standard  input 


CHARACTER  Dir*80 , Header*20 , Index*80 
- Read  the  directory  name  header. 
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non  on  on  oooooo  on  on  non  n  n  on 


READ  (*, 5,ERR=30,END=20)  Header 
5  FORMAT  (A20) 

10  FORMAT  (A80) 

IF  (Header (1: 6) .NE. 'LOOKUP' )  GOTO  30 
READ  (*, 10,ERR=30,END=20)  Dir 

- Read  the  index  file  header. 

READ  (*, 5,ERR=30,END=20)  Header 

IF  (Header (1: 17) .NE. 'LOOKUP  INDEX  FILE')  GOTO  30 

READ  (*, 10,ERR=30,END=20)  Index 

RETURN 

- Report  some  error  and  bail  out. 

0  Header  =  'EOF' 

0  WRITE  (*,40)  Dir , Index, Header 

0  FORMAT ( 'ERROR :  expecting  a  header,  followed  by  a  string:'/ 

>  5x, 'DIRECTORY  NAME'/  A60/  5x, 'INDEX  FILE  NAME'/  A60/ 

>  5x, 'with : ' , A20 ) 

STOP 

END 

SUBROUTINE  Rdrh (Rhhi , Rh) 


Read  the  relative  humidity  measurements  into  an  array. 


CHARACTER  Header *20 
INTEGER  N 

REAL  Rhhi (100) ,Rh(100) 

- Skip  rh  title  line. 

READ  (*, 10,ERR=50,END=40)  Header 
10  FORMAT  (A20) 

IF  ( Header (1: 16 ) .NE. 'HUMIDITY  PROFILE')  GOTO  50 

- --Read  one  humidity  datum. 

N  =  1 

0  READ  (*, *,ERR=30,END=40)  Rhhi (N) , Rh (N) 

- Negative  height  flag  means  end  of  the  list;  returns. 

IF  (Rhhi(N) .LT.0.0)  RETURN 

- Check  that  the  humidity  is  between  1  and  99  percent. 

IF  (Rh(N) .LE.l.O.OR.Rh(N) .GT.99.9)  GOTO  30 

_ Ready  to  fill  next  element,  then  over-write  last  until  neg 

IF  (N.LT.lOO)  N  =  N+1 
GOTO  20 

- Error  for  unreasonable  humidity . 

30  WRITE  (*,*)  'ERROR:  Rh  must  be  in  the  range  1  to  99  percent.' 
WRITE  (*,*)  'Instead,  f ound : ' , Rh (N) , '  N= ' , N 
STOP 

- Other  kinds  of  missing  or  bad  input  for  humidity. 

40  Header  =  'EOF' 

50  WRITE  (*,*)  'ERROR:  expected  humidity  prof ile . ' 

WRITE  (*,*)  'Instead,  f ound :', Header 

STOP 

END 

SUBROUTINE  Rdtc ( Tchi , Tc ) 


Read  the  temperature  measurements  into  an  array. 


C 


CHARACTER  Header *20 
INTEGER  N 

REAL  Tchi(lOO) ,Tc(100) 
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C  - Temperature  profile  title  line. 

READ  ( * , 10 , ERR=50 , END=40 )  Header 
10  FORMAT  (A20) 

IF  (Header (1:19) .NE. 'TEMPERATURE  PROFILE')  GOTO  50 
C 

C  - Read  one  temperature  datum. 

N  =  1 

20  READ  (*, *,ERR=30,END=40)  Tchi (N) , Tc (N) 

C 

0  - Negative  height  means  end  of  the  list. 

IF  (Tchi (N) .LT. 0 . 0)  RETURN 
C 

C  - Check  that  the  temperature  is  reasonable. 

IF  (Tc(N) .LT. -50.0. OR. Tc(N) .GT. 50.0)  GOTO  30 


_ Ready  to  fill  next  element.  Over-write  #100  until  neg .  height. 

IF  (N.LT.lOO)  N  =  N+1 
GOTO  20 


C 

0  - Error  for  unreasonable  temperature. 

30  WRITE  (*,*)  'ERROR:  Tc  must  be  in  the  range  -50  C  to  50  C. 
WRITE  (*,*)  'Instead,  f ound : ' , Tc (N) , '  N= ' , N 
STOP 


0  _ Other  kinds  of  missing  or  bad  input  for  temperature. 

40  Header  =  'EOF' 

50  WRITE  (*,*)  'ERROR:  expected  temperature  profile. 

WRITE  {*,*)  'Instead,  f ound :', Header 

STOP 

END 

SUBROUTINE  Rdtob(Tob) 

- - 

C  Read  a  one-third  octave  band  spectrum,  bands  0-43. 

- - 

INTEGER  I 
REAL  Tob(0:43) 

C 

C  - Read  the  spectrum,  bands  0  -  ,9,  10  -  19,  etc. 

READ  (*, *,ERR=100)  (Tob ( I ) , 1=0 , 43 ) 

RETURN 

C 

100  WRITE  {*,*)  'ERROR:  bad  number  in  spectrum' 

STOP 

END 

SUBROUTINE  Rdver 

- - 

C  Read  the  version  number  for  Lookup. 

- - 

CHARACTER  Header *20 
REAL  Ver 


0  - Read  the  title  line  and  check  if  we've  got  the  right  input. 

READ  (*, 10,ERR=40,END=30)  Header 
10  FORMAT  (A20) 

IF  (Header (1:19) .NE. 'LOOKUP  LEVELS  INPUT')  GOTO  40 
C 

C  - Check  the  version  number. 

Ver  =0.0 

READ  ( * , * , ERR=2  0 , END=3  0 )  Ver 
C 

C  - Return  if  we  have  the  right  input  file. 

IF  (Ver.EQ.1.0)  RETURN 
C 

0  - Error  for  bad  file  version. 

20  WRITE  (*,*)  'ERROR:  expecting  Lookup  program  version:  1.0.' 
WRITE  (*,*)  'Instead,  found:', Ver 
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STOP 

C 

- Error  for  other  bad  or  missing  input. 

3  0  Header  =  'EOF'  ■.  j  / 

40  WRITE  (*,*)  'ERROR:  expecting  Lookup  input  header. 

WRITE  (*,*)  'Instead,  f ound Header 

STOP 

END 

SUBROUTINE  Rdwswd {Whi , Wx, Wy ) 

C  Read  the  wind  speed  and  direction  measurements  into  _ 

CHARACTER  Header *20 

INTEGER  N  , ,  , 

REAL  Pi , Whi ( 100 ) , Wdir ( 100 ) , Ws ( 100 ) , Wx ( 10 0 ) , Wy ( 10 0 ) 

C 

C  - Compute  Pi. 

Pi  =  4 . 0*ATAN(1 . 0) 

C 

Q  - Skip  a  the  wind  profile  title  line. 

READ  (*, 10,ERR=50,END=40)  Header 
10  FORMAT  (A20) 

IF  (Header (1: 12) .NE. 'WIND  PROFILE')  GOTO  50 

Q  - Read  one  wind  speed,  wind  direction  datum. 

N  =  1 

20  READ  ( * , * , ERR=3 0 , END=40 )  Whi (N) , Ws (N) , Wdir (N) 

C 

Q  - Negative  height  means  end  of  the  list;  returns. 

IF  (Whi(N) .LT.0.0)  RETURN 

C  _ Don't  allow  negative  wind  speeds  or  speeds  above  40  m/s. 

IF  (Ws(N) .LT. 0.0. OR. Ws(N) .GT. 40.0)  GOTO  30 
C 

Q  - Use  elements  1  to  99,  only. 

IF  (N.LT.lOO)  THEN 

Q  - Convert  to  Wx,  and  Wy.  (Wdir  is  meteo.  wind  direction!!) 

Wx(N)  =  -Ws (N) *SIN(Wdir (N) /ISO . 0*Pi) 

Wy(N)  =  -Ws (N) *COS (Wdir (N) /ISO . 0*Pi) 

C 

Q  - Ready  to  fill  next  element. 

N  =  N+1 
END  IF 
GOTO  2  0 
c 

U  - Error  for  unreasonable  wind  speed  or  illegal  numbers. 

30  WRITE  (*,*)  'ERROR:  Ws  must  be  in  the  range  0  to  40  m/s' 

WRITE  (*,*)  'Instead,  found:  Ws (N) = ' , Ws (N) , '  N= ' , N 
STOP 

Q  - Other  kinds  of  missing  or  bad  input  for  wind. 

40  Header  =  'EOF' 

50  WRITE  {*,*)  'ERROR:  expected  wind  profile.' 

WRITE  (*,*)  'Instead,  f ound :', Header 

STOP 

END 
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APPENDIX  B:  Sample  Input  File 


LOOKUP  LEVELS  INPUT  FILE,  VERSION: 

1.0 

LOOKUP  DIRECTORY  NAME: 

/usr/ white/ Prog /Lookup/ 

LOOKUP  INDEX  FILE: 
indexlO.fil 
SOURCE  HEIGHT  (m) 

11.7 

RECEIVER  HEIGHT  (m) 

2.1 

RECEIVER  RANGE  (m)  (0=ALL) 

0.0 

RECEIVER  AZIMUTH  (deg) 

90.0 

GROUND  CLASSIFICATION  (HARD=0,  POROUS=l) 

0.51 

GROUND  ROUGHNESS  HEIGHT  (m) 

0  1 

TEMPERATURE  PROFILE;  HEIGHT (m)  TEMPERATURE  (deg  C) 

0.1  9.2 

1.010.0 

1.4  10.4 

20.0  5.3 
100.03.3 
200.0  -1.0 
-1.0  -1.0  THE  END 

HUMIDITY  PROFILE;  HEIGHT  (m)  REL .  HUMIDITY  (%) 

0.5  50.0 

1.0  55.0 

10.0  45.0 
100.0  55.0 

200.0  20.0 

-1.0  -1.0  THE  END  ,  , 

WIND  PROFILE;  HEIGHT  (m) ,  WIND  SPEED  (m/s),  DIRECTION  (deg) 

2.0  4.2  240.0 

5.0  6.0  300.0 

-1.0  -1.0  -1.0  THE  END 

SOURCE  LEVEL  SPECTRUM  (dB) ;  1/3  OCTAVE  BANDS  0-43 


0.0 

0  . 

0 

0  . 

0 

0  . 

0 

0  . 

0 

70.0 

73 

.0 

76 

.0 

79 

.0 

82 

.0 

85.0 

88 

.0 

91 

.0 

94 

.0 

97 

.0 

100.0 

97 

.0 

94 

.0 

91 

.0 

88 

85.0 

82 

.0 

79 

.0 

76 

.0 

73 

.0 

70.0 

0  . 

0 

0  . 

0 

0  . 

0 

0  . 

0 

0.0 

0  . 

0 

0  . 

0 

0  . 

0 

0  . 

0 

0.0 

0  . 

0 

0  . 

0 

0  . 

0 

0. 

0 

0.0 

0  . 

0 

0  . 

0 

0  . 

0 

REFERENCE  DIST  (m) 

250.0  „ 

USER  WEIGHTING  SPECTRUM  (dB) ;  1/3  OCTAVE  BANDS  0  43 


0  . 

0 

0  . 

0 

0  . 

0 

0  . 

0 

0  . 

0 

0  . 

0 

0  . 

0 

0  . 

0 

0  . 

0 

0  . 

0 

10 

.0 

10 

.0 

10 

.0 

10 

.0 

10 

.0 

0  . 

0 

0  . 

0 

0. 

0 

0  . 

0 

0  . 

0 

0. 

0 

0  . 

0 

0  . 

0 

0  . 

0 

0  . 

0 

0  . 

0 

0. 

0 

0. 

0 

0  . 

0 

0. 

0 

0. 

0 

0  . 

0 

0. 

0 

0  . 

0 

0. 

0 

0. 

0 

0  . 

0 

0  . 

0 

0  . 

0 

0  . 

0 

0. 

0 

0  . 

0 

0. 

0 

0  . 

0 
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LOOKUP  LEVELS  INPUT  FILE,  VERSION: 

1.0 

LOOKUP  DIRECTORY  NAME: 

/usr /white /Prog /Lookup/ 

LOOKUP  INDEX  FILE: 
indexlO . f il 
SOURCE  HEIGHT  (m) 

11.7  - 

.RECEIVER  HEIGHT  (m) 

2.1 

RECEIVER  RANGE  (m)  (0=ALL) 

10.0 

RECEIVER  AZIMUTH  (deg) 

90 . 0 

GROUND  CLASSIFICATION  (HARD=0,  POROUS=l) 

0.51 

GROUND  ROUGHNESS  HEIGHT  (m) 

0.1  , 

TEMPERATURE  PROFILE;  HEIGHT (m)  TEMPERATURE  (deg  C) 

0.1  9.2 

1.0  10.0 

1.4  10.4 

20.0  5.3 
100.0  3.3 

200.0  -1.0 

-1.0  -1.0  THE  END 

HUMIDITY  PROFILE;  HEIGHT  (m)  REL .  HUMIDITY  (%) 

0.5  50.0 

1.0  55.0 

10.0  45.0 
100.0  55.0 

200.0  20.0 

-1.0  -1.0  THE  END  ,  ^  ^ 

WIND  PROFILE;  HEIGHT  (m) ,  WIND  SPEED  (m/s) ,  DIRECTION  (deg) 

2.0  4.2  240.0 

5.0  6.0  300.0 

-1.0  -1.0  -1.0  THE  END 

SOURCE  LEVEL  SPECTRUM  (dB) ;  1/3  OCTAVE  BANDS  0-43 


0.0 

0  . 

0 

0. 

0 

0  . 

0  , 

0  . 

0 

70.0 

73 

.0 

76 

.0 

79 

.  0 

82 

.0 

85.0 

88 

.  0 

91 

.0 

94 

.  0 

97 

.0 

100 . 0 

97 

.0 

94 

.0 

91 

.0 

88 

85.0 

82 

.0 

79 

.0 

76 

.0 

73 

.0 

70.0 

0  . 

0 

0  . 

0 

0  . 

0 

0  . 

0 

0.0 

0  . 

0 

0  . 

0 

0  . 

0 

0  . 

0 

0  .  0 

0  . 

0 

0  . 

0 

0  . 

0 

0  . 

0 

0.0 

0  . 

0 

0  . 

0 

0  . 

0 

REFERENCE  DIST  (m) 

250.0 

USER  WEIGHTING  SPECTRUM  (dB) ;  1/3  OCTAVE  BANDS  0-43 


0  . 

0 

0  . 

0 

0  . 

0 

0  . 

0 

0  . 

0 

0  . 

0 

0  . 

0 

0  . 

0 

0  . 

0 

0  . 

0 

10 

.0 

10 

.0 

10 

.0 

10 

.0 

10 

.0 

0. 

0 

0  . 

0 

0  . 

0 

0  . 

0 

0  . 

0 

0  . 

0 

0  . 

0 

0  . 

0 

0  . 

0 

0  . 

0 

0  . 

0 

0  . 

0 

0  . 

0 

0  . 

0 

0  . 

0 

0. 

0 

0  . 

0 

0  . 

0 

0  . 

0 

0  . 

0 

0  . 

0 

0  . 

0 

0. 

0 

0  . 

0 

0  . 

0 

0  . 

0 

0  . 

0 

0. 

0 

0  . 

0 
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APPENDIX  C:  Sample  Index  File 


ahnnl4' 

0 

-14 

bhnnl4' 

0 

-14 

chnnl4 ' 

0 

-14 

dhnnl4' 

0 

-14 

asnnl4 ' 

0 

-14 

bsnnl4 ' 

0 

-14 

csnnl4 ' 

0 

-14 

dsnnl4 ' 

0 

-14 

ahnnl2 ' 

0 

-12 

bhnnl2 ' 

0 

-12 

chnnl2 ' 

0 

-12 

'  dhnnl2 ' 

0 

-12 

'  asnnl2 ' 

0 

-12 

’  bsnnl2 ' 

0 

-12 

'  csnnl2 ' 

0 

-12 

'  dsnnl2 ' 

0 

-12 

'  ahnnlO ' 

0 

-10 

'bhnnlO ' 

0 

-10 

'  chnnlO ' 

0 

-10 

'  dhnnlO ' 

0 

-10 

'  asnnlO ' 

0 

-10 

'bsnnlO ' 

0 

-10 

'  csnnlO ' 

0 

-10 

'  dsnnlO ' 

0 

-10 

'  ahnnOS ' 

0 

-8 

'bhnnOS ' 

0 

-8 

'  chnnOS ' 

0 

-8 

'dhnnOS ' 

0 

-8 

'  asnnOS ' 

0 

-8 

'bsnnOS ' 

0 

-8 

'  csnnOS ' 

0 

-8 

'dsnnOS' 

0 

-8 

'ahnn06 ' 

0 

-6 

'bhnn06 ' 

0 

-6 

'chnn06' 

0 

-6 

'dhnn06 ' 

0 

-6 

' asnnO  6 ' 

0 

-6 

'bsnn06 ' 

0 

-6 

' csnnOS ' 

0 

-6 

' dsnn06 ' 

0 

-6 

' ahnnOB ' 

0 

-5 

'bhnnOB ' 

0 

-5 

' chnnOB ' 

0 

-5 

' dhnnOB ' 

0 

-5 

' asnnOB ' 

0 

-5 

'bsnnOB ' 

0 

-5 

' csnnOB ' 

0 

-5 

' dsnnOB ' 

0 

-5 

' ahnn04 ' 

0 

-4 

'bhnn04 ' 

0 

-4 

' chnn04 ' 

0 

-4 

'dhnn04 ' 

0 

-4 

' asnn04 ' 

0 

-4 

'bsnn04 ' 

0 

-4 

' csnn04 ' 

0 

-4 

' dsnn04 ' 

0 

-4 

' ahnn02 ' 

0 

-2 

'bhnn02 ' 

0 

-2 

' chnn02 ' 

0 

-2 

'dhnn02' 

0 

-2 

' asnn02 ' 

0 

-2 

'bsnn02 ' 

0 

-2 

0 

2 

1.2 

0 

5 

1.2 

0 

10 

1.2 

0 

100 

1.2 

1 

2 

1.2 

1 

5 

1.2 

1 

10 

1.2 

1 

100 

1.2 

0 

2 

1.2 

0 

5 

1.2 

0 

10 

1.2 

0 

100 

1.2 

1 

2 

1.2 

1 

5 

1.2 

1 

10 

1.2 

1 

100 

1.2 

0 

2 

1.2 

0 

5 

1.2 

0 

10 

1.2 

0 

100 

1.2 

1 

2 

1.2 

1 

5 

1.2 

1 

10 

1.2 

1 

100 

1.2 

0 

2 

1.2 

0 

5 

1.2 

0 

10 

1.2 

0 

100 

1.2 

1 

2 

1.2 

1 

5 

1.2 

1 

10 

1.2 

1 

100 

1.2 

0 

2 

1.2 

0 

5 

1.2 

0 

10 

1.2 

0 

100 

1.2 

1 

2 

1.2 

1 

5 

1.2 

1 

10 

1.2 

1 

100 

1.2 

0 

2 

1.2 

0 

5 

1.2 

0 

10 

1.2 

0 

100 

1.2 

1 

2 

1.2 

1 

5 

1.2 

1 

10 

1.2 

1 

100 

1.2 

0 

2 

1.2 

0 

5 

1.2 

0 

10 

1.2 

0 

100 

1.2 

1 

2 

1.2 

1 

5 

1.2 

1 

10 

1.2 

1 

100 

1.2 

0 

2 

1.2 

0 

5 

1.2 

0 

10 

1.2 

0 

100 

1.2 

1 

2 

1.2 

1 

5 

1.2 
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'csnn02  ' 

0 

-2 

1 

10 

1.2 

'  dsnn02 ' 

0 

-2 

1 

100 

1.2 

'  ahnnOl ' 

0 

-1 

0 

2 

1.2 

'  bhnnOl ' 

0 

-1 

0 

5 

1.2 

'  chnnOl ' 

0 

-1 

0 

10 

1.2 

'  dhnnOl ' 

0 

-1 

0 

100 

1.2 

'asnnOl' 

0 

-1 

1 

2 

1.2 

'  bsnnOl ' 

0 

-1 

1 

5 

1.2 

' csnnOl ' 

0 

-1 

1 

10 

1.2 

' dsnnO 1 ' 

0 

-1 

1 

100 

1.2 

'ahnpOO ' 

0 

0 

0 

2 

1.2 

'bhnpOO ' 

0 

0 

0 

5 

1 . 2 

'  chnpO  0 ' 

0 

0 

0 

10 

1 . 2 

'  dhnpOO ' 

0 

0 

0 

100 

1.2 

' asnpOO ' 

0 

0 

1 

2 

1 . 2 

'bsnpOO' 

0 

0 

1 

5 

1 . 2 

'  csnpO  0 ' 

0 

0 

1 

■  10 

1 . 2 

'dsnpOO ' 

0 

0 

1 

100 

1 .2 

'  ahnpOl ' 

0 

1 

0 

2 

1 . 2 

'bhnpOl ' 

0 

1 

0 

5 

1 . 2 

' chnpOl ' 

0 

1 

0 

10 

1 . 2 

'dhnpOl ' 

0 

1 

0 

100 

1 . 2 

'asnpOl ' 

0 

1 

1 

2 

1 . 2 

'bsnpOl ' 

0 

1 

1 

5 

1 . 2 

' csnpO 1 ' 

0 

1 

1 

10 

1 . 2 

'dsnpOl ' 

0 

1 

1 

100 

1.2 

' ahnp02 ' 

0 

2 

0 

2 

1 . 2 

' bhnp02  ' 

0 

2 

0 

5 

1 . 2 

' chnp02 ' 

0 

2 

0 

10 

1 . 2 

' dhnp02  ' 

0 

2 

0 

100 

1 . 2 

' asnp02 ' 

0 

2 

•  1 

2 

1 . 2 

'bsnp02 ' 

0 

2 

1 

5 

1 . 2 

' csnp02 ' 

0 

2 

1 

10 

1 . 2 

' dsnp02  ' 

0 

2 

1 

100 

1 . 2 

' ahnp04 ' 

0 

4 

0 

2 

1 . 2 

'bhnp04 ' 

0 

4 

0 

5 

1 . 2 

' chnp04 ' 

0 

4 

0 

10 

1 . 2 

' dhnp04 ' 

0 

4 

0 

.  100 

1.2 

'  asnpO'4 ' 

0 

4 

1 

2 

1 . 2 

'bsnp04 ' 

0 

4 

1 

5 

1.2 

' csnp04 ' 

0 

4 

1 

10 

1.2 

' dsnp04 ' 

0 

4 

1 

100 

1.2 

' ahnpO  5 ' 

0 

5 

0 

2 

1.2 

'bhnp05' 

0 

5 

0 

5 

1 . 2 

' chnp05 ' 

0 

5 

0 

10 

1.2 

'dhnpOB' 

0 

5 

0 

100 

1.2 

' asnp05 ' 

0 

5 

1 

2 

1 . 2 

' bsnpO  5  ' 

0 

5 

1 

5 

1 . 2 

' csnpO  5 ' 

0 

5 

1 

10 

1 . 2 

' dsnpOB ' 

0 

5 

1 

100 

1.2 

' ahnp06 ' 

0 

6 

0 

2 

1 . 2 

'bhnp06 ' 

0 

6 

0 

5 

1 . 2 

' chnp06 ' 

0 

6 

0 

10 

1 . 2 

'dhnp06 ' 

0 

6 

0 

100 

1.2 

' asnp06 ' 

0 

6 

1 

2 

1 . 2 

'bsnp06 ' 

0 

6 

1 

5 

1.2 

' csnp06 ' 

0 

6 

1 

10 

1.2 

' dsnpOS ' 

0 

6 

1 

100 

1.2 

' ahnpOS ' 

0 

8 

0 

2 

1 . 2 

'bhnpOS' 

0 

8 

0 

5 

1.2 

' chnpOS ' 

0 

8 

0 

10 

1 . 2 

' dhnpOS ' 

0 

8 

0 

100 

1.2 

' asnpOS ' 

0 

8 

1 

2 

1 . 2 

'bsnpOS ' 

0 

8 

1 

5 

1 . 2 

' csnpOS ' 

0 

8 

1 

10 

1 . 2 

' dsnpOS ' 

0 

8 

1 

100 

1.2 
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'  ahnplO ' 

0 

10 

'bhnplO ' 

0 

10 

' chnplO  ' 

0 

10 

' dhnplO  ' 

.  0 

10 

' asnplO  ' 

0 

10 

'bsnplO  ' 

0 

10 

' csnplO  ' 

0 

10 

'dsnplO  ' 

0 

10 

' ahnpl2  ' 

0 

12 

'bhnpl2 ' 

0 

12 

'chnpl2' 

0 

12 

' dhnpl2 ' 

0 

12 

' asnpl2 ' 

0 

12 

'bsnpl2 ' 

0 

12 

' csnpl2 ' 

0 

12 

' dsnpl2 ' 

0 

12 

' ahnpl4 ' 

0 

14 

'bhnpl4' 

0 

14 

' chnpl4 ' 

0 

14 

' dhnpl4 ' 

0 

14 

' asnpl4 ' 

0 

14 

' bsnpl4 ' 

0 

14 

' csnpl4 ' 

0 

14 

' dsnpl4 ' 

0 

14 

' ahgnl4 ' 

1 

-14 

' bhgnl4 ' 

1 

-14 

' chgnl4 ' 

1 

-14 

' dhgnl4 ' 

1 

-14 

' asgnl4 ' 

1 

-14 

'bsgnl4 ' 

1 

-14 

' csgnl4 ' 

1 

-14 

' dsgnl4 ' 

1 

-14 

' ahgnl2 ' 

1 

-12 

'bhgnl2' 

1 

-12 

' chgnl2 ' 

1 

-12 

'dhgnl2 ' 

1 

-12 

'asgnl2 ' 

1 

-12 

'bsgnl2 ' 

1 

-12 

' csgnl2 ' 

1 

-12 

' dsgnl2 ' 

1 

-12 

' ahgnlO ' 

1 

-10 

'bhgnlO ' 

1 

-10 

' chgnlO ' 

1 

-10 

'dhgnlO' 

1 

-10 

' asgnlO ' 

1 

-10 

'bsgnlO ' 

1 

-10 

' csgnlO ' 

1 

-10 

' dsgnlO ' 

1 

-10 

' ahgnOS ' 

1 

-8 

'bhgnOS' 

1 

-8 

'chgnOS' 

1 

-8 

'dhgnOS' 

1 

-8 

' asgnOS ' 

1 

-8 

'bsgnOS' 

1 

-8 

'csgnOS' 

1 

-8 

' dsgnOS ' 

1 

-8 

'ahgn06 ' 

1 

-6 

'bhgn06 ' 

1 

-6 

' chgn06 ' 

1 

-6 

'dhgn06 ' 

1 

-6 

' asgn06 ' 

1 

-6 

'bsgnOS ' 

1 

-6 

' csgnO  6 ' 

1 

-6 

' dsgn06 ' 

1 

-6 

'ahgnOB ' 

1 

-5 

'bhgnOB ' 

1 

-B 

0 

2 

1.2 

0 

5 

1.2 

0 

10 

1.2 

0 

100 

1.2 

1 

2 

1.2 

1 

B 

1.2 

1 

10 

1.2 

1 

100 

1.2 

0 

2 

1.2 

0 

B 

1.2 

0 

10 

1.2 

0 

100 

1.2 

1 

2 

1.2 

1 

B 

1.2 

1 

10 

1.2 

1 

100 

1.2 

0 

2 

1.2 

0 

B 

1.2 

0 

10 

1.2 

0 

100 

1.2 

1 

2 

1.2 

1 

B 

1.2 

1 

10 

1.2 

1 

100 

1.2 

0 

2 

1.2 

0 

B 

1.2 

0 

10 

1.2 

0 

100 

1.2 

1 

2 

1.2 

1 

5 

1.2 

1 

10 

1.2 

1 

100 

1.2 

0 

2 

1.2 

0 

5 

1.2 

0 

10 

1.2 

0 

100 

1.2 

1 

2 

1.2 

1 

5 

1.2 

1 

10 

1.2 

1 

100 

1.2 

0 

2 

1.2 

0 

B 

1.2 

0 

10 

1.2 

0 

100 

1.2 

1 

2 

1.2 

1 

5 

1.2 

1 

10 

1.2 

1 

100 

1.2 

0 

2 

1.2 

0 

B 

1.2 

0 

10 

1.2 

0 

100 

1.2 

1 

2 

1.2 

1 

B 

1.2 

1 

10 

1.2 

1 

100 

1.2 

0 

2 

1.2 

0 

5 

1.2 

0 

10 

1.2 

0 

100 

1.2 

1 

2 

1.2 

1 

B 

1.2 

1 

10 

1.2 

1 

100 

1.2 

0 

2 

1.2 

0 

B 

1.2 

33 


'  chgnO  5 ' 

1 

-5 

0 

10 

1.2 

'  dhgnOB ' 

1 

-5 

0 

100 

1.2 

'  asgnOS ' 

1 

-5 

1 

2 

1 . 2 

'  bsgnO  5 ' 

1 

-5 

1 

5 

1 . 2 

' csgnOB ' 

1 

-5 

1 

10 

1 . 2 

' dsgnO  5 ' 

1 

-5 

1 

100 

1 . 2 

' ahgn04 ' 

1 

-4 

0 

2 

1 . 2 

' bhgn04 ' 

1 

-4 

0 

5 

1 . 2 

' chgn04 ' 

1 

-4 

0 

10 

1 . 2 

' dhgn04 ' 

1 

-4 

0 

100 

1 . 2 

' asgn04 ' 

1 

-4 

1 

2 

1 . 2 

'bsgn04 ' 

1 

-4 

1 

5 

1 . 2 

' csgn04 ' 

1 

-4 

1 

10 

1 . 2 

'dsgn04'  ' 

1 

-4 

1 

100 

1 . 2 

' ahgn02 ' 

1 

-2 

0 

2 

1 . 2 

'bhgn02 ' 

1 

-2 

0 

5 

1 .  z 

' chgn02 ' 

1 

-2 

0 

10 

1 . 2 

' dhgn02 ' 

1 

-2 

0 

100 

1 . 2 

1  o 

' asgn02 ' 

1 

-2 

1 

2 

1  .  z 

I  n 

'bsgn02 ' 

1 

-2 

1 

5 

1 .  z 

"I  o 

' csgn02 ' 

1 

-2 

1 

10 

1 . 2 

' dsgn02 ' 

1 

-2 

1 

100 

1 . 2 

1  o 

'ahgnOl ' 

1 

-1 

0 

2 

1  .  Z 

1  o 

'bhgnOl ' 

1 

-1 

0 

5 

1  .  z 

' chgnO 1 ' 

1 

-1 

0 

10 

1 . 2 

'dhgnOl' 

1 

-1 

0 

100 

1 . 2 

' asgnOl ' 

1 

-1 

1 

2 

1  .  Z 

' bsgnO 1 ' 

1 

-1 

1 

5 

X  •  z 

' csgnOl ' 

1 

-1 

1 

10 

1.2 

'dsgnOl ' 

1 

-1 

1 

100 

1 . 2 

' ahgpOO ' 

1 

0 

0 

2 

1 .  z 

1  o 

'bhgpOO ' 

1 

0 

0 

5 

1 . 2 

' chgpOO ' 

1 

0 

0 

10 

1 .2 

'dhgpOO ' 

1 

0 

0 

100 

1.2 

'asgpOO ' 

1 

0 

1 

2 

1 . 2 

'bsgpOO ' 

1 

0 

1 

5 

1 , 2 

' csgpOO ' 

1 

0 

1 

10 

1 . 2 

'dsgpOO ' 

1 

0 

1 

100 

1 . 2 

'ahgpOl ' 

1 

1 

0 

2 

1 . 2 

'bhgpOl ' 

1 

1 

0 

5 

1 .2 

'chgpOl' 

1 

1 

0 

10 

1 . 2 

' dhgpOl ' 

1 

1 

0 

100 

1.2 

' asgpOl ' 

1 

1 

1 

2 

1 . 2 

'bsgpOl ' 

1 

1 

1 

5 

1 .2 

' csgpOl ' 

1 

1 

1 

10 

1 . 2 

' dsgpOl ' 

1 

1 

1 

100 

1.2 

' ahgp02 ' 

1 

2 

0 

2 

1 . 2 

'bhgp02 ' 

1 

2 

0 

5 

1 . 2 

' chgp02 ' 

1 

2 

0 

10 

1 . 2 

'dhgp02 ' 

1 

2 

0 

100 

1 . 2 

' asgp02 ' 

1 

2 

1 

2 

1 .2 

'bsgp02 ' 

1 

2 

1 

5 

1 .2 

' srcc/csgp02 .n ' 

1 

2 

1 

10  1.2 

'dsgp02 ' 

1 

2 

1 

100 

1 .2 

'ahgp04 ' 

1 

4 

0 

2 

1 . 2 

'bhgp04 ' 

1 

4 

0 

5 

1 . 2 

' chgp04 ' 

1 

4 

0 

10 

1 . 2 

'dhgp04 ' 

1 

4 

0 

100 

1 .2 

'asgp04 ' 

1 

4 

1 

2 

1 . 2 

' bsgp04 ' 

1 

4 

1 

5 

1 . 2 

' csgp04 ' 

1 

4 

1 

10 

1 .2 

' dsgp04 ' 

1 

4 

1 

100 

1 .2 

' ahgpOB ' 

1 

5 

0 

2 

1 . 2 

'bhgpOB ' 

1 

5 

0 

5 

1 .2 

' chgpOS ' 

1 

5 

0 

10 

1 .2 

' dhgpOB ' 

1 

5 

0 

100 

1.2 
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'asgpOB ' 

1 

5 

1 

2 

1.2 

'bsgpOB ' 

1 

5 

1 

5 

1 . 2 

'  csgpOS  '■ 

1 

5 

1 

10 

1 . 2 

'dsgpOB' 

1 

5 

1 

100 

1 . 2 

1  o 

' ahgp06 ' 

1 

6 

0 

2 

1 . 2 
-1  o 

'bhgp06 ' 

1 

6 

0 

5 

1 . 2 

' chgp06 ' 

1 

6 

0 

10 

1 . 2 

'dhgpOG' 

1 

6 

0 

100 

1 .2 

1  o 

' asgp06 ' 

1 

6 

1 

2 

1 . 2 

1  o 

'bsgp06' 

1 

6 

1 

5 

1 . 2 

' csgpO  6 ' 

1 

6 

1 

10 

1 .2 

' dsgp06 ' 

1 

6 

1 

100 

1 . 2 

1  0 

' ahgpOS ' 

1 

8 

0 

2 

1 . 2 

1  o 

'bhgpOS' 

1 

8 

0 

5 

1 . 2 

1  o 

' chgpOS ' 

1 

8 

0 

10 

1 . 2 

'dhgpOS ' 

1 

8 

0 

100 

1 . 2 

1  o 

' asgpOS ' 

1 

8 

1 

2 

1 . 2 

1  O 

'bsgpOS ' 

1 

8 

■  1 

5 

1 . 2 

' csgpOS ' 

1 

8 

1 

10 

1 . 2 

' dsgpOS ' 

1 

8 

1 

100 

1 . 2 

1  0 

'ahgplO ' 

1 

10 

0 . 

2 

J.  .  2 

1  O 

'bhgplO ' 

1 

10 

0 

5 

1.2 

' chgplO ' 

1 

10 

0 

10 

1 . 2 

' dhgplO ' 

1 

10 

0 

100 

1 . 2 

1  0 

' asgplO ' 

1 

10 

1 

2 

X  •  2 

1  o 

'bsgplO ' 

1 

10 

1 

5 

i  .  2 

I  o 

'  csgplO ' 

1 

10 

1 

10 

1 . 2 

' dsgplO ' 

1 

10 

1 

100 

1 . 2 

1  0 

'ahgpl2 ' 

1 

12 

0 

2 

X  •  2 

'bhgpl2 ' 

1 

12 

0 

5 

X  .  2 

I  o 

' chgpl2 ' 

1 

12 

0 

10 

1 . 2 
•1 

'dhgpl2 ' 

1 

12 

0 

100 

1 . 2 

1  0 

' asgpl2 ' 

1 

12 

1 

2 

X  .  2 

1  o 

'bsgpl2' 

1 

12 

1 

5 

1 . 2 
■1  o 

'csgpl2 ' 

1 

12 

1 

10 

1 . 2 

' dsgpl2 ' 

1 

12 

1 

100 

1 . 2 

'ahgpl4 ' 

1 

14 

0 

2 

X  .  2 
1  o 

'bhgpl4 ' 

1 

14 

0 

5 

1 . 2 

' chgpl4 ' 

1 

14 

0 

10 

1 . 2 

' dhgpl4 ' 

1 

14 

0 

100 

1 . 2 

'asgpl4 ' 

1 

14 

1 

2 

1 . 2 
*1  o 

'bsgpl4 ' 

1 

14 

1 

5 

1 . 2 

'  csgpl4 ' 

1 

14 

1 

10 

1 . 2 

'dsgpl4 ' 

1 

14 

1 

100 

1 . 2 
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APPENDIX  D:  Sample  Output  File 


BEGIN  RANGE,  FLAT 

A, 

C 

AND 

USER 

1 

.  0 

133 

.  0 

103  . 

,3 

130 

.0 

139  , 

.0 

1 

.3 

133 

.  0 

103  . 

,3 

130 

.  0 

138  , 

.9 

1 

.  6 

133 

.0 

103  . 

,2 

130 

.0 

138  , 

,9 

2 

.  0 

132 

.9 

103  . 

,2 

129 

.9 

138, 

.8 

2 

.  5 

132 

.8 

103  . 

,  1 

129 

.8 

138, 

.7 

3 

.2 

132 

.6 

102  . 

,9 

129 

.6 

138  . 

.  6 

4 

.0 

132 

.4 

102  . 

,7 

129 

.4 

138  . 

,3 

5 

.0 

132 

.0 

102  . 

,3 

129 

.0 

138  . 

.  0 

6 

.3 

131 

.5 

101. 

,8 

128 

.5 

137. 

,  5 

7 

.9 

130 

.8 

101. 

.1 

127 

.8 

136. 

.8 

10 

.  0 

129 

.8 

100  . 

.  1 

126 

.8 

135, 

.7 

12 

.  6 

128 

.7 

98  . 

.9 

125 

.  6 

134, 

.  6 

15 

.8 

127 

.3 

97  . 

.  5 

124 

.3 

133  , 

.2 

20 

.0 

125 

.7 

96  . 

.0 

122 

.7 

131 , 

.7 

25 

.  1 

123 

.9 

94. 

.2 

120 

.9 

129  , 

.9 

31 

.  6 

122 

.2 

92  . 

.4 

119 

.1 

128  , 

.  1 

39 

.8 

120 

.3 

90  , 

.6 

117 

.3 

126, 

.2 

50 

.  1 

118 

.3 

88. 

.5 

115 

.3 

124, 

.2 

63 

.  1 

116 

.3 

86, 

.  6 

113 

.3 

122  , 

.3 

79 

.4 

114 

.3 

84  , 

.5 

111 

.3 

120  , 

.2 

100 

.  0 

112 

.2 

82  , 

.5 

109 

.2 

118., 

.1 

125 

.9 

110 

.1 

80  , 

.4 

107 

.1 

116, 

.0 

158 

.5 

107 

.9 

78, 

.2 

104 

.9 

113  , 

.8 

199 

.5 

105 

.7 

76. 

.0 

102 

.7 

111, 

.7 

251 

.2 

103 

.4 

73  , 

.7 

100 

.4 

109, 

.4 

316 

.2 

101 

.  1 

71 , 

.4 

98 

.1 

107  , 

.  1 

398 

.1 

98 

.7 

69, 

.0 

95 

.7 

104, 

.7 

501 

.2 

96 

.2 

66, 

.5 

93 

.2 

102, 

.2 

631 

.0 

93 

.5 

63  , 

.8 

90 

.5 

99  , 

.5 

794 

.3 

90 

.7 

61, 

.0 

87 

.7 

96, 

.7 

1000 

.0 

87 

.7 

58, 

,0 

84 

.7 

93  , 

.7 

1258 

.9 

84 

.4 

54, 

,7 

81 

.4 

90  , 

.4 

1584 

.9 

80 

.8 

51 . 

,1 

77 

.8 

86, 

.8 

1995 

.3 

76 

.7 

47  . 

,0 

73 

.7 

82, 

.7 

2511 

.9 

72 

.1 

42  , 

,4 

69 

.1 

78  , 

.1 

3162 

.3 

6  6 

.9 

37  , 

.2 

63 

.9 

72, 

.9 

3981 

.1 

60 

.8 

31, 

.1 

57 

.8 

66, 

.8 

5011 

.9 

53 

.6 

23  . 

.9 

50 

.6 

59  , 

.  6 

6309 

.6 

51 

.4 

19, 

.6 

48 

.3 

57  , 

.5 

7943 

.3 

49 

.3 

17  , 

.0 

46 

.2 

55, 

.5 

10000 

.  0 

47 

.2 

14, 

.5 

44 

.0 

53  , 

.4 

12589 

.3 

45 

.1 

11, 

.8 

41 

.9 

51 

.3 

15848 

.9 

43 

.  0 

,  9 , 

.1 

39 

.7 

49 

.3 

19952 

.  6 

40 

.8 

6, 

.4 

37 

.4 

47 

.2 

36 


BEGIN  RANGE,  FLAT,  A,  C  AND  USER 

10.0  129.8  100.1  126.8  135.7 


25. 

0 

25  . 

0 

25  . 

0 

25  . 

0 

95  . 

0 

98  . 

,0 

101 . 

0 

104  . 

,0 

110  . 

,0 

113  . 

,0 

116. 

,0 

119  , 

.0 

125. 

.0 

122  . 

.  0 

119  . 

.0 

116. 

.0 

110  . 

,  0 

107  , 

.0 

104. 

,0 

101, 

.0 

95, 

.0 

25, 

.0 

25, 

.0 

25, 

.0 

25, 

.0 

25 

.0 

25, 

.0 

25 

.0 

25 

.  0 

25 

.0 

25 

.0 

25 

.0 

25 

.0 

25 

.0 

25 

.0 

25 

.0 

25 . 0 
107.0 
122.0 
113.0 
98.0 
25.0 
25.0 
25.0 
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ECA  Team  Distribution 


Chief  of  Engineers 
ATTN;  CEHEC-IM-LP  (2) 

ATTN:  CEHEC-IM-LH  (2) 

ATTN;  CERD-L 
ATTN:  CEMP-CE 
ATTN:  CEMP-EA 
ATTN:  CEMP-EI  (2) 

ATTN;  CEMP-ZA 
ATTN:  CEMP-ZM  (2) 

ATTN;  DAIM-EO 

HQ  USAF/LEEEU  20332 

US  Army  Europe 
ODCS/Engineer  09014 
ATTN:  AEAEN-FE 
ATTN:  AEAEN-ODCS 

AMC  22333 
ATTN:  AMCEN-A 

Fort  Belvoir,  VA  22060 
ATTN:  Water  Resource  Center 
ATTN;  CECC-R 
ATTN:  NACEC-FB 

Picatinny  Arsenal  07801 
ATTN:  Library 

US  Military  Academy  10996 
ATTN:  Facilities  Engineer 
ATTN:  Dept  of  Geography  & 
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